REVIEW ARTICLE 



Dynamics of vortices and interfaces 
in superfluid "^He 

A. P. Finne^, V. B. Eltsov^ 2, R. Hannineni'", 
N. B. Kopnin^ J. Kopu\ M. Krusius\ 
M. Tsubota'', and G. E. Voloviki ^ 

^Low Temperature Laboratory, Helsinki University of Technology, P.O. Box 
2200, FIN-02015 HUT, Finland 

^Kapitza Institute for Physical Problems, 119334, Moscow, Russia 
^Landau Institute for Theoretical Physics, 119334, Moscow, Russia 
''Department of Physics, Osaka City University, Osaka 558-8585, Japan 



Abstract. Rapid new developments have occurred in superfluid hydrodynamics 
since the discovery of a host of unusual phenomena which arise from the diverse 
structure and dynamics of quantized vortices in ^He superfluids. These have 
been studied in rotating flow with NMR measurements which at best provide an 
accurate mapping of the different types of topological defects in the superfluid 
order parameter fleld. Four observations are reviewed here: (1) the interplay 
of different vortex structures at the first order interface between the two major 
superfluid '^He phases, '^He-A and '^He-B; (2) the shear flow instability of 
this phase boundary, which is now known as the superfluid Kelvin-Helmholtz 
instability; (3) the hydrodynamic transition from turbulent to regular vortex 
dynamics as a function of increasing dissipation in vortex motion; and (4) 
the peculiar propagation of vortex lines in a long rotating column which even 
in the turbulent regime occurs in the form of a helically twisted vortex state 
behind a well-developed vortex front. The consequences and implications of these 
observations are discussed, as inferred from measurements, numerical calculations, 
and analytical work. 
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1. Introduction 

1.1. Helium superfluids 

Until 1972 the only known example of a truly inviscid fluid was superfluid *He-II at 
low flow velocities. Today its primacy is challenged by the discovery of superfluid 
^He in 1972 and the gaseous Bose-Einstein condensates in 1995. Nevertheless, in low 
temperature physics ^He-II remains the epitome of a superfluid, the benchmark to 
which to compare to. Its perfect inviscid flow is known to persist only at velocities 
below some low critical limit at which quantized vortex lines are formed. These 
are topologically stable linear defects in the superfluid order parameter field with 
extraordinary properties. One of them is their turbulent flow, customarily known 
as superfluid turbulence or quantum turbulence which appears, for instance, 
when the applied flow velocity is suddenly increased well above the critical limit. 
Recently the dynamics of quantized vortex lines has gained renewed interest, activated 
by the differences and similarities which have been discovered while exploring the 
other superfluid systems, both superfluid '^He and Bose-Einstein condensates. A 
recent review by Vinen and Niemela summarizes with updated references our 
understanding on vortex dynamics and turbulence in ^He-II. 

Here we are concerned with the superfluid phases of '^He which usher into 
superfluid hydrodynamics a broad spectrum of new phenomena. These are associated 
with the structure and dynamics of quantized vortices and other topological defects. 
As bulk liquid, superfluid '^He can exist in three different phases, of which the 
two major phases are ^He-A and ^He-B. The third phase, ''He-Ai, exists at high 
magnetic fields around the zero- field superfluid transition temperature Tc . The flow 
properties of '^He-B are isotropic in the absence of external magnetic fields, resembling 
those of ■*He-II with its quantized vorticity. In contrast, ^He-A is highly anisotropic 
and the most extraordinary superfluid of all that we know. Its applications as a 
model system in physics have far reaching implications 0]. This review deals with 
recent observations on vortex dynamics, primarily in '^He-B, which have been made 
in uniformly rotating flow with noninvasive NMR measurement. Related reviews can 
be found in Refs. 0017]. 

Although '^He-II and '^He-B both are isotropic helium superfluids and in many 
respects rather similar, if compared to ^He-A, nevertheless, important differences 
prevail which have profound influence on the resulting superfluid hydrodynamics. The 
implications from these differences have been appreciated only lately. One of them 
concerns the vortex-core radius, whose length scale in both cases is determined by 
the coherence length £,{T) of the superfluid state. In *He-II the core radius is of 
atomic scale ^ 0.1 nm, while in ■^He-B it is > 10 nm and thus at least two orders 
of magnitude larger. This difference is not simply quantitative but has substantial 
impact on the interactions of the vortex with the container wall, on critical velocities, 
vortex formation, and surface pinning. The implications from this difference became 
apparent in the first half of the 1990ies, when singic-vortcx formation was observed 
in large open volume measurements in "^He-B, but only in flow through micron-size 
orifices in '^He-II . 

The second major difference is the viscosity of the normal component in the two 
superfluids. In •^Hc-B it has oil-like viscosity and is practically always in a state of 
laminar flow. In contrast, the normal component of ^He-II is one of the least viscous 
fluids known. Its flow becomes easily turbulent, which in turn influences the flow of the 
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superfluid component, resulting in complicated mutual- friction coupled turbulence of 
the normal and superfluid fractions. Thus the absence of turbulence in the flow of the 
normal component of '^He-B amounts to a considerable simplification at temperatures 
above the zero temperature limit (where the normal component is present). 

The third difference which influences profoundly the dynamics of vortices is 
mutual friction dissipation, the damping which takes place when a vortex moves with 
respect to the normal component. In the Fermi superfluid '^He-B mutual friction 
between the vortex and the normal component is mediated by fermionic quasiparticle 
states in the vortex core |H| El 1101 the so-called fermion zero modes. Their 
properties are described by a theory similar to the BCS theory of superconductivity, 
according to which the dimensionless temperature-dependent parameter q(T), which 
is the ratio of the dissipative and reactive components of the mutual friction force, 
is a function of the normalized temperature T/Tc and depends almost exponentially 
on temperature. It crosses unity at around T ~ 0.6 Tc. At temperatures above 
this division point the number of vortices is generally found to remain constant 
in dynamic processes. In contrast, at lower temperatures vortices become easily 
unstable in externally applied flow which causes an increase in the vortex number 
owing to superfluid turbulence. In comparison, in the Bose liquid ^He-II mutual 
friction dissipation is small in the usual regime of measurements and vortex dynamics 
is practically always turbulent. The regular vortex number conserving flow might be 
expected only within a few /xK below the superfluid transition temperature T\. From 
this temperature regime there are no experiments available on vortex dynamics yet. 
Even there, the low viscosity of the normal component might cause the coupled flow 
to become turbulent. 

In addition to their hydrodynamic differences, '^He-II and ■^He-B experiments 
often use different techniques to create and detect vortex lines. The temperature 
required for superfluid '^He is a factor of 10^ lower than for ^He-II. This sets restrictions 
on the type of experiments that can be conducted on '^He superfluids. Uniform 
rotation can be used in any temperature range to create counterflow of the normal 
and superfluid components. Owing to better control over vortex formation in "^He 
superfluids, rotation has there proven to be a useful means to apply flow. 

As for vortex detection, in '^He superfluids nuclear magnetic resonance (NMR) 
measurement provides a practical noninvasive method to count the number of vortex 
lines and to study their dynamics. In '^He-B NMR methods can be used from 
down to about 0.2 Tc, while in ^He-A measurements at even lower temperatures 
should be possible. As the temperature decreases measurements based on equilibrium 
state techniques become increasingly less sensitive. This is the case also in ^He-B 
NMR, where changes in the order parameter texture from superfluid counterflow, 
vortices, and other control parameters gradually vanish or saturate with decreasing 
temperature. This is an unfortunate constraint, since today the T ^ limit is of 
great interest in superfluid hydrodynamics. 

In the zero-temperature limit, where the normal component becomes exponent- 
ially rarefied, the only measuring methods developed so far for the study of vortices 
in ^He-B employ vibrating wires ^^j) spheres ^Hli or grids jj^. These resonantly 
oscillating objects can be employed as sensitive sensors of their hydrodynamic 
environment in a quiescent He bath, for instance to create and detect vortices. The 
oscillation is driven at amplitudes where the flow velocity at the surface of the vibrating 
body exceeds the critical value for Cooper-pair breaking |15) . In the zero-temperature 
regime of ballistic quasiparticle motion, a second resonant sensor oscillating at low 



Dynamics of vortices and interfaces in superBuid ^He 



5 



drive in the linear regime can then be used to track deviations in the exponentiaUy 
temperature-dependent equilibrium quasiparticle density or the quasiparticles 
scattered from the flow field around a vortex or a tangle of vortices These 
techniques have turned out to provide efficient new tools for vortex studies ^Sl and 
are now in the forefront of future research. They can also be used for constructing 
ultra-sensitive dark matter detectors ^H] ■ The lack of suitable measuring techniques 
has also been an obstacle in vortex studies of *He-II at the lowest temperatures. A 
promising new development is here the use of micron-size charged vortex rings for the 
analysis of different vortex states. With this method both rotating arrays of rectilinear 
lines and turbulent tangles can be distinguished and monitored |20|. 

1.2. Novel phenomena in superfluid ^ He hydrodynamics 

Here we outline briefly the four main topics which are the subject of this review. 
The flrst is concerned with the fundamental difference in the structure of quantized 
vortex lines in the A and B phases of superfluid '^He. This set of questions is peculiar to 
superfluid ■^He. It is the only presently known system where vortices can be studied at a 
stable flrst order interface between two coherent states which belong to the same order 
parameter manifold. Here the phase of the order parameter is continuous across the 
interface and thus vortices can, in principle, cross the interface continuously. This is 
quite unlike other interfaces, for instance between phase-separated layers of superfluid 
^He above a solution of '^He in superfluid ''He. In this latter case the quantized vortices 
in the two layers belong to different superfluid systems and can end at the interface 
with little relation to each other. 

^He-A is an anisotropic liquid where, in a typical experimental situation in a 
magnetic fleld, the vortex core is formed on a length scale which is at least three orders 
of magnitude larger than in ^He-B. This scale is not set by the pairing interaction, but 
by the tiny dipolar coupling between the spin and orbital momenta of the Cooper pairs. 
The structural length scale of quantized vorticity is not the superfluid coherence length 
£,{T) > 10 nm, but the healing length ^d{T) > 10 /im associated with the dipolar spin- 
orbit coupling. The typical A-phase vortex is doubly quantized, i.e. its circulation is 
twice that of the "^He-B vortex. This difference between the vortices poses a problem 
when they interact at the interface between these two superfluids in a rotating sample: 
How is the large core doubly-quantized A-phase vortex matched with the small core 
singly-quantized B-phase vortex across the AB interface? Measurements elucidating 
this question led to the unexpected observation of dissipationless shear flow between 
the two superfluids at the AB interface. The stability issue of this superfluid shear-flow 
state is one the topics discussed in this review. 

The possibility of constructing the shear flow state arises from the different 
conditions of vortex formation in the two superfluids, owing to the large difference 
in vortex core radius. The core of the ^He-B vortex is intermediate between that 
in "^He-II and '^He-A, which leads to important consequences. On one hand, being 
larger than the microscopic core in ^He-II, pinning and surface roughness at bounding 
walls is not as important as in '^He-II. With carefully chosen and prepared container 
surfaces pinning sites can be avoided, so that pinned remnant vortices do not exist. In 
such cases, substantial vortex-free flow can be reached in a cylindrical rotating sample, 
before intrinsic vortex formation starts to intervene at relatively high critical velocities. 
In contrast, in ^He-II vortcx-free flow has generally little practical meaning because, 
even at very low velocities, remnant vorticity leads to efficient vortex formation. An 
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important exception is flow through a sub-micron-size aperture in a thin membrane 
where vortices are swept away from the immediate vicinity of the high-velocity flow 
and do not have a chance to become pinned there |21j . 

On the other hand, the core radius of the '^He-B vortex is much smaller than 
that of a continuous vortex in ■^He-A. As a result, the critical velocity for intrinsic 
vortex formation in ^He-B is much larger than in '^He-A. This makes it possible 
to prepare a flow state in which vortices are already forming on the A-phase side 
of the AB interface, while on the B-phase side the vortex-free irrotational Landau 
state persists. Such a situation leads to a shear-flow state in which the superfluid 
components of the two superfluids are sliding with respect to each other at the AB 
interface. The relative flow of the two superfluids is frictionless and, for the first 
time, provides a perfect arrangement for the experimental investigation of the classical 
Kelvin-Hclmholtz instability which was theoretically predicted hundred and fifty years 
ago (Sec. 12. 4|) . The reason for this unique situation is that in conventional viscous 
liquids the threshold for the Kelvin-Helmholtz instability, where the formation of 
surface waves or ripplons on the interface starts, is always obscured by the influence 
of viscosity. 

Unexpectedly, even in the perfect superfluid conditions, the critical velocity of the 
AB-interface instability does not match the classical result derived for ideal inviscid 
fluids (Sec. |3Jl. However, a modified criterion for the onset of the instability proved 
to be in excellent agreement with the superfluid experiments, although it appears 
to lead to paradoxical consequences at first glance. This instability threshold is not 
determined by the velocity Vs2— Vsi of the relative superfluid motions in the two liquids 
across the interface, but the instability would occur even if the two liquids would have 
the same velocity or if there is a single superfluid with a free surface. These new 
features result from the two-fluid nature of the superfluid liquid, from the presence 
of the superfluid and normal fractions. The instability threshold is determined by 
the velocities Vsi_2 — Vn of each superfluid with respect to the reference frame of 
the container walls and thus with respect to the normal fractions of the two liquids, 
which in thermodynamic equilibrium move together with the walls. The free surface 
of a superfluid bath with respect to its gas phase (or vacuum at the low temperatures) 
becomes unstable, when in the reference frame of the normal component, the superfluid 
velocity reaches the critical threshold value |22|2S1- In the case of several superfluid 
fractions («) in the same liquid, such as neutron and proton superfluids in a neutron 
star, the threshold is determined by some combination of the superfluid velocities 

Surprisingly, the superfluid Kelvin-Helmholtz instability has many features in 
common with the instability of quantum vacuum beyond the event horizon or even in 
the ergoregion of the black hole. The ergoregion is defined as the region at the interface 
where the energy of surface waves, or ripplons, is negative. The ripplon excitations of 
the AB interface also provide a connection to the presently popular idea in cosmology, 
according to which matter in our Universe is confined to hypersurfaces, which are 
multidimensional membranes, or branes, in a multidimensional space. Branes can be 
represented by topological defects, such as domain walls and strings, and by interfaces 
between different quantum vacua. In our case, the brane is defined by the AB interface 
between two quantum vacua - the two superfluid phases of '^He. The instability of the 
AB interface is in one-to-one correspondence to the instability of quantum vacuum in 
the brane world. It occurs in the ergoregion because of the interaction between the 
matter on the brane (represented by ripplons) and the matter in higher-dimensional 
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space (represented by quasiparticles in bulk superfluids). 

Measurements on the AB interface instability also revealed new properties about 
superfluid turbulence (SecQ)), which is the third main topic of this review. The later 
nonlinear stage of the AB interface instability results in the injection of a tight bundle 
of small vortex loops in the rapidly flowing vortex-free ^He-B. It was found that the 
injection leads to turbulence in "^He-B at temperatures below a critical onset value 
Ton ~ 0.6 Tc. The temperature of this hydrodynamic transition turned out to depend 
only on the dimensionless intrinsic parameter q{T), the ratio of the dissipative and 
reactive mutual friction components. The fact that in superfluid hydrodynamics a 
transition to turbulence occurs as a function of mutual friction dissipation at g ~ 1 
was discovered for the first time in ^He-B j^SI- It divides the evolution of the injected 
vortices to regular vortex number conserving dynamics at q > 1 and to turbulence 
^■t 9 ^ 1- The main reason why the transition has not been observed in superfluid 
^He and only in '^He-B is the favorable range of values of the parameter q{T). This 
fortunate coincidence, the presence of the hydrodynamic transition in the middle of 
the experimentally accessible temperature range, makes it possible to explore the 
dynamics in both flow regimes under otherwise similar conditions. 

Of particular significance has been the exponentially steep temperature depend- 
ence of the mutual friction dissipation. It has allowed a whole new genre of studies 
on how turbulence switches on, when one or a few vortices which are far apart, are 
introduced in vortex-free flow. A further unexpected phenomenon is the evolution and 
propagation of the vortices in a long rotating cylinder or column after the turbulence 
has switched on. It turns out that the propagation takes the form of a spiralling 
vortex front which travels longitudinally and rotates azimuthally with respect to the 
cylinder walls and thereby expands in the unstable vortex-free state. Behind the front 
an ordered helically twisted vortex bundle forms where the vortices are in a force-free 
configuration. This twisted state is already close in energy to the final state of solid- 
body rotation, to which it relaxes when the vortex front has reached the end plate of 
the rotating cylinder (SecE}. Thus the front separates here in effect the metastable 
vortex-free Landau state from the equilibrium vortex state. The motion of the front 
and the helically twisted state can be monitored with the NMR measurement. These 
observations and their interpretation form the fourth topic of the present review. They 
are not a special characteristic of superfluid ■^He-B, but apply equally to superfluid 
''He-II, for instance. However, they became possible in ^He-B because of better control 
over vortex formation and the possibility to create vortex-free flow at relatively high 
flow rate. A further characteristic of these measurements is a longer sample cylinder 
than has been used before in rotating measurements, with two separate detectors, 
which made it possible to record changes in the flow state in different parts of the 
rotating column as a function of time. 

The above issues have been in the forefront of recent research and are in the focus 
of this review. They demonstrate new features of superfluid hydrodynamics and often 
arise owing to the multi-component order parameter of the '^He superfluids. In such 
cases they cannot be reproduced with the 'classical' ''He-II superfluid. However, some 
of these phenomena or their analogues might be present in the new superfluid states 
of gaseous bosonic or fermionic atom clouds. 
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2. Hydrodynamics of rotating helium superfluids 

2.1. Helium superfluids and vortex lines 

Customarily ''He-II is described with, a wave function = e*''', where (p is the 
phase factor. The superfluid velocity is then defined as the gradient of the phase, 
Vg — (K/27r)V(/3 where k = 2Trh/m4. Since the curl of a gradient vanishes identically, 
V X Vg = 0, the bulk superflow is irrotational. In principle rotational flow is thus 
excluded, but by forming quantized vortex lines the condensate can accommodate to 
rotating flow. In its simplest form a line vortex is a stable string-like object with 
a central hard core where the order parameter vanishes in the center, thus forming 
a line singularity in the coherent order-parameter field. Since the condensate phase 
changes by 2?™ around the core, where n is an integer number, the circulation of the 
superfluid velocity around the vortex core is quantized: 

(j) dr ■ -Vs — UK , (1) 

and K plays the role of the circulation quantum. This persistent superfluid current 
around the core stores kinetic energy, providing the vortex with an energy per unit 
length, or line tension, which equals 

e.^lps £ dr{2nr)v! = ^ In (^) • (2) 

Here the upper (rv) and lower (re) cutoffs are determined by the inter- vortex distance 
and the core size of the order of the superfluid coherence length ^, respectively. Vortex 
lines with a singly-quantized structure n — I are thus energetically favourable. The 
superfluid hydrodynamics which follows from the introduction of the quantized vortex 
lines has been described in textbooks |26l I27| . 

The order parameter in superfluid ^He relates to the wave function of the Cooper 
pairs, and has a complicated internal structure. Nevertheless, in ■^He-B it still contains 
an explicit phase variable ip. The above considerations remain valid with the exception 
that the circulation quantum is now given by k = 21:^1/ {2ra^) — 0.066 mm^/s, where 
2m3 is the mass of the two '^He atoms in the Cooper pair, rather than the single atom 
mass 7714 in the '*He-II circulation quantum k — 0.099 mm^/s. In what follows we 
use the same symbol k to denote the circulation quantum in both "^He and ^He, with 
appropriate values for each particular fluid. The superfluid order parameter does not 
vanish in the '^He-B vortex core, but the order parameter state in the core is different 
from that in the bulk. Two different vortex-core structures are known to exist: an 
axisymmetric core at high temperatures and high pressures, and a nonaxisymmctric 
at low temperatures |28| |29j . A first-order phase transition, which under equilibrium 
conditions occurs at 0.60 at 29 bar pressure |2HI) separates these two core structures. 
At low pressures ^' ^ 15 bar only the nonaxisymmctric core exists. In both cases the 
core radius Tc is approximately equal to the coherence length ^ ^ 10 nm. An interesting 
curiosity to note is that this transition was the first phase transition ever observed 
within a defect, when it was discovered in 1981 30 . A third vortex structure in bulk 
^He-B is the spin-mass vortex, a combination of a linear and a planar defect with both 
spin and mass flow currents around its core |^ I32L |^ . It will not be discussed in this 
review. 
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Figure 1. (Left) Macroscopic structure of quantized vortex line in He superfluids. 
The core radius Tc is on the order of the superfluid coherence length in *He-II 
~ 0.1 nm) and ^He-B > 10 nm), but in '^He-A the length scale is the healing 
length of the dipolar spin-orbit interaction (^jj > 10 /im). (Right) Orbital order 
parameter texture of the soft core of the double-quantum vortex in ^He-A in 
magnetic field. The cones indicate the local direction and rotation of the orbital 
order parameter triad of unit vectors 1, ih, n. The topological winding number 
of the 1 texture is n = 2. The texture is nonaxisymmetric: it is composed of a 
circular half, or meron, and a hyperbolic meron, each with 2it circulation. 



^He-II and "^He-B are traditional Landau superfluids in that their superflow is 
potential, V x Vg = 0, unless vortex line defects are present. In •^He-A, where the 
phase and the orbital structure (represented by the orbital vector 1) of the order 
parameter are linked together, this condition is no longer strictly satisfied. Instead, 
the so-called Mermin-Ho relation holds : 

V X Vs -^^'^ijkJiiyij X Vffe). (3) 

This implies that rotational superfluid flow can be accomplished via an inhomogeneous 
order-parameter texture l(r). However, the energy cost of the necessary spatial 
variations, resulting from the rigidity of the order parameter, gives rise to a finite 
critical superflow velocity also in this system. At this velocity, vorticity with 
continuously winding structure of the order parameter orientation is formed so that 
in most cases no hard vortex core is involved. In the simplest form the structure of an 
isolated continuous vortex has the following spatial distribution of the orbital 1-field: 

l(p, 0) = z cos ri{p) + p sin ■q{p) . (4) 

Here z, p and (f) are the unit vectors of the cylindrical coordinate system; 77 (p) changes 
from ry(0) = to »y(oo) = tt. This winding 1 texture forms the so-called continuous 
soft core of the vortex j34| . since it is in this region where the non-zero vorticity of 
superfluid velocity is concentrated: 

Vs(/5,0) = - cos77(p)] ^ , Vxvs = -^sin?7 (|^) z .(5) 

The circulation of the superfluid velocity around a contour enclosing the soft-core 
region is quantized, f dr ■ Vs = 2k, corresponding to the quantization number n = 2. 
Thus the object described by Eq. (0} is a continuous double-quantum vortex. By 
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following the 1 field across the cross section of the soft-core texture, it is noted that 
the 1 vector goes through all possible orientations on the unit sphere. Such a topology 
of the vortex cross section in two spatial dimensions is known as a skyrmion. 

In the magnetic field of the NMR measurements the continuous vortex is deformed 
and its structure is nonaxisymmetric, see Fig. ^ However, its topology is robust to 
deformations, and the circulation remains the same: ^ dr • Vg = 2k. It is important 
to note that, since even in the soft-core region the order parameter retains its bulk 
structure, the core size of the continuous ■^He-A vortex is not set by the coherence 
length ^ 10 nm of the superfluid state. Instead, the relevant length scale is the three 
orders of magnitude larger dipolar healing length ^ 10 /im which originates from 
the spin-orbit coupling. 

Using the two halves of the skyrmion texture, the circular and hyperbolic merons 
(Fig.^ as basic building blocks, other structures of continuous vorticity can be formed. 
An example are the various continuous periodic vortex textures in zero or low magnetic 
field Another important structure is the vortex sheet [23 ES] which competes 

for living space with the double-quantum vortex line. A concise lexicon of these various 
structures can be found in Ref. pSj . 

The concept of the quantized vortex line dates back to Lars Onsager (1949) 

and Richard Feynman (1955) @J who found that the Landau irrotationality 
requirement V x Vg = has to be lifted at singular lines where V x Vg ^ 0. In the 
case of ^He-A these principles were put to a severe test which they finally survived 
when, in the context of the work of Mermin and Ho in 1976 ^HIj Chechetkin (1976) 
02] and Anderson and Toulouse (1977) ^IH^ the first example of a continuous vortex 
texture was proposed. 

2.2. Vortex states in rotating superfluid 

The identification of the vortex structures of superfluid ■^He, and of the phase 
transitions separating these different structures, is based to a large extent on NMR 
measurements on a rotating sample. In rotation the vorticity V x Vg is aligned parallel 
to the rotation axis Jl and generally forms a regular array over the cross section of the 
cylindrical sample. This is a particularly simple situation where both the structural 
and dynamic properties of these vortex structures can be analyzed. 

The minimum energy configuration in rotation is the state with the equilibrium 
number of rectilinear vortex lines A^cq, which on average mimics solid-body rotation 
of the superfluid, i.e. (vg) = Jl x r, or (V x Vg) = 217. Since (V x Vg) = nKu^, 
the vortex density in the bulk is — 2Q/{nK,). The formation of a new vortex 
is associated with an energy barrier that has to be overcome before an elementary 
vortex loop can be nucleated. At sufficiently low applied flow velocities this is not 
possible, and metastable states with a vortex number TV smaller than iVoq can be 
formed. These consist of a central vortex cluster (Fig. EJ, with any number of 
vortex lines < A < A'oq. Within the cluster the rectilinear lines are packed to 
their equilibrium density = 2fl/(nK), confined by the counterfiow of the normal 
and superfiuid components which circulates around the cluster with the velocity 
V = Vn — Vg = [ri/' — ukN/ (27rr)] cf>. The first term is the velocity of the normal 
component, locked to co-rotation with the cylindrical container (with radius R), while 
the second term arises from the combined persistent superfiow of the N rectilinear 
vortex lines in the central cluster. An extreme case is the Landau state - the vortex- 
free state with A^ = and Wg = (as expressed in the rest frame of the laboratory). 
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Figure 2. Schematic representation of a vortex cluster confined by the 
azimuthally circulating counterflow of the normal and superfluid components to 
the center of the rotating sample. The areal density of rectihnear singly-quantized 
vortex lines within the cluster is 2Q / k and thus their number in a cluster of radius 
Rc is N = ttRI 2C/k. 



This is the state of maximum kinetic energy in the rotating frame. In many of the 
rotating experiments described below it is the initial state, the starting point for 
the measurements. Independently of N ^ the maximum counterflow velocity is at the 
cylindrical wall at r = i?. This we call the velocity of the externally applied flow or 
the rotation drive of the cylindrical rotating container. 

At constant rotation the stationary states are thus the equilibrium vortex state 
and the various metastable states with a depleted vortex cluster. In an ideal cylinder, 
which is exactly aligned parallel to the rotation axis, it is possible to have more 
than the equilibrium number of vortices A^oq, owing to a finite annihilation barrier 
01] . Experimentally the exact value of iVcq is important for calibrating the measuring 
signals from a state with a well-defined configuration and number of vortices. Transient 
time-dependent rotating states are created in accelerating or decelerating rotation 
|45|. In Sees. 0] - we describe measurements where rotation is kept constant and 
the dynamics evolves from vortex seed loops which have been introduced by external 
means into initially vortex-free counterflow. 

2. 3. Critical velocity of vortex formation 

The lowest critical velocity in a rotating superfluid is that at which the free energy of a 
rectilinear vortex line flrst becomes negative in the container frame. The corresponding 
angular velocity is known as the Feynman critical velocity, fici — K/(27ri?^) \n{R/rc) 
|1T] . It is analogous to the critical held Hci for type II superconductors. For a rotating 
cylinder with a radius of a few mm, f2ci ~ 10~^ rad/s and is thus very small. However, 
although at O > Od it becomes energetically favorable to introduce vortices in the 
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Rotation velocity Q, (rad/s) 



Figure 3. Measurement of single vortex formation as a function of the applied 
rotation velocity f2 at high temperatures in ^He-B 1471 . (Top) The vertical axis 
shows the NMR absorption in the Larmor region of the '^He-B spectrum. Vortex 
formation starts with the first step-like increase, but the critical threshold at Qc 
is identified from the third step where the critical flow velocity in the bottom 
plot reaches a more stable value. (Bottom) Counterfiow velocity at the cylinder 
wall V = QR — kN/{2itR), where N is the number of steps already formed. The 
maximum possible value of counterflow for this sample container is defined by the 
horizontal dashed line, v = = QcR- 



rotating sample, some mechanism for their formation is required. A number of such 
mechanisms exist, owing to sources both extrinsic and intrinsic to the superfluid itself, 
each with its characteristic critical velocity. For more discussion see Ref. 46_. Here 
we only summarize the basic ideas important for the overall picture. 

In practice, in ^He-B the lowest critical velocity fie, which controls the formation 
of vortices, is found to depend on the shape and size of the container and the roughness 
of its surfaces. The simplest and most ideal case is a smooth- walled cylinder which 
is mounted with its symmetry axis as parallel to the rotation axis as possible. The 
surface quality is dependent on the choice of material, the fabrication of the seams in 
the corners, and the cleanliness of the walls. Even residual gases, such as air or water, 
will condense on the wall during cool down, form small crystallites, and may determine 
ric- In a good sample cylinder of typical radius 2 ~ 3 mm, 0,^. is relatively high, of 
order 1 - 4 rad/s, so that large vortex-free counterflow can be achieved before the first 
vortex is formed. In the worst case extrinsic sources govern Q^- For instance, it can 
be determined by some pinning site, a piece of dirt, at which a remnant vortex may 
remain pinned indefinitely. If this site is occupied and is increased to the critical 
value associated with the site, the pinned remnant loop will start to evolve. 

In the most favorable case flc arises from a combination of extrinsic and intrinsic 
reasons, if vortex formation takes place at a sharp surface asperity in the form of 
a pointed spike 03 1321 ■ At a very sharp spike the local velocity can exceed the 
average velocity at the wall by one to two orders in magnitude. Thus superfluidity 
will be broken first at this location when is increased to fie, and a small vortex 
loop is formed 0H! • The loop then evolves to a rectilinear vortex line and reduces the 
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counterflow velocity at the cylinder wall to a sub-critical value v = flcR — k/{2ttR). 
If Q is increased further by external means, vortex formation occurs recurrently at 
the same site every time when the counterflow reaches the critical value Vc — flcR- 
Here Vc is therefore the limit for vortex-free flow in this container. An example of 
such vortex formation in single-quantum events as a function of f2 is illustrated by the 
staircase pattern in Fig. 13 This measurement has been performed in a container with 
R = 2 mm |47| . A similar measurement with *He-II has been demonstrated only with 
flow through orifices of sub- micron size |^ 011 ■ 

An estimate of the intrinsic critical velocity, and of the energy barrier which 
inhibits the formation of an elementary vortex loop, can be obtained from the following 
simple consideration ^H]. The barrier is determined by the smallest possible vortex 
ring. Since the radius of such a ring cannot be smaller than the core radius (of the 
order of ^), the energy of the smallest ring can be estimated from Eq. (jSJ as £' ~ Psn"^^- 
This gives -E/^bT ~ {^/ a){T-p /T), where we have used Ps ~ m/a^ for the superfluid 
density and T-p = fi^ /2ma?k-Q ~ IK for the degeneracy temperature of the quantum 
Fermi liquid, with a as the interatomic distance. For '^He-B we obtain E/k-oT > 10^. 
This should be compared to a similar estimate E/k-oT > 1 for "'He-II, where the core 
size and the coherence length are £, ^ a. 

How to overcome such an energy barrier j5()ll51| ? The rate for thermal activation 
over the barrier is oc expi— E / k^T) , and thus the barrier is practically impenetrable 
by thermal activation or quantum tunneling |52| at the appropriate temperatures for 
^He superfluids (T ~ 1 mK) 0^1 QZ] • Both mechanisms become effective only when 
the local velocity at the asperity reaches a value extremely close to the threshold where 
the energy barrier disappears, and the hydrodynamic instability of the flow occurs. 
This occurs at a velocity of order Vc ^ k/^. In ^He-B, this critical velocity for vortex 
formation is comparable to the Landau critical velocity for quasiparticle creation - 
the pair-breaking velocity Vph = A/pp ~ k/^, where A is the superfluid energy gap. 

In contrast, in ^He-A the smallest possible vortex loop is of the order of the 
soft-core radius, the healing length of the spin-orbital coupling ^ 10 /xm. This 
is several orders of magnitude larger than the coherence length ^. Consequently, the 
critical velocity for A-phase vortex formation, Vc ^ k/£.d, is considerably smaller |46) 
while the energy barrier is higher than in the B phase. 

Therefore, in practical experimental conditions neither thermal activation nor 
quantum tunneling are of importance in ''He superfluids. Instead, vortex formation 
takes place when the average counterflow velocity at the wall is increased to the point 
where the local velocity at the sharpest asperity reaches the critical value, the barrier 
vanishes, and the process thus becomes an instability. In principle, pair breaking and 
quasiparticle emission might occur already at a slightly lower velocity than when 
the barrier actually disappears, and this might finally trigger the hydrodynamic 
instability, which then results in vortex formation. The process might happen in 
the following manner: near the asperity the local velocity reaches the pair breaking 
value, the creation and emission of quasiparticles increases the density of the normal 
component, and as a result ps decreases. Due to the conservation of current pgVg the 
superfluid velocity then increases, enhancing the radiation of quasiparticles, which 
increases Vg further. The final result from the development of such a hydrodynamic 
instability will be vortex formation. However, whatever is the real mechanism of 
the instability generated by the flow in the vicinity of a protuberance, it limits the 
maximum counterflow velocity that can be achieved in a given sample container. With 
careful preparation of the surfaces, the critical velocity Vc has been raised up to about 
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0.1 - 0.4 Wpb in cylinders from fused quartz with i? ^ 3 mm. 

The situation is quite different in ''He-II. Ahhough the maximum possible 
superfluid velocity, the Landau value, is three orders of magnitude higher, the 
nucleation barrier height is 1 - 10 K and comparable to the ambient temperature 
of IK. In the flow through sub-micron-size orifices thermal activation has been 
found to be an important mechanism in vortex nucleation [531 149j (at the lowest 
temperatures even quantum tunneling has been argued to exist jK31l55j V In contrast, 
in applications of bulk volume *He-II flow it is assumed that there always exist an 
abundance of remnant vortex loops pinned to walls |56| which start to expand in low 
applied flow. 

In '^He-B surface pinning is expected to be much less important than in ^He-II 
because the vortex core radius is more than two orders of magnitude larger. In the 
best conditions in a clean quartz cylinder there are no pinned remnant vortices. In 
such cases no other kind of information about surface pinning exists at present time. 
However, dynamic remnant vortices are present at low temperatures. When mutual 
friction dissipation becomes exponentially small, it may take hours for the last vortex 
to annihilate at the container wall in conditions of zero applied flow. If a remnant 
vortex loop happens to be around when flow is reapplied, it starts to evolve and may 
generate any number of new independent vortex loops |57| . In fact, this situation 
is expected to prevail also in bulk ''^He-II over most of the experimentally accessible 
temperature range and thus no source with abundant pinned remnant vortices is 
needed to create large numbers of vortices and turbulent flow. 

In '^He-B we thus expect that a genuine intrinsic critical velocity is determined 
by the most effective instability, since the vortex formation barrier is impenetrable at 
all temperatures and velocities below an upper limit Vc{T,P). This feature has been 
utilized to study the different instabilities described in the later sections. The upper 
limit is a container specific critical velocity (which may change from one cool down 
to another, depending on the container's preparation), of which an example is seen in 
Fig. 13 In such a measurement the criterion for vortex formation is the lowest critical 
velocity, in other words it is the vortex structure with the lowest Vc which is formed. 
If on the other hand one wants to establish the true equilibrium vortex structure, 
one has to slowly cool the sample at constant flow velocity below Tc- At Tc critical 
velocities vanish and the criterion for the selection becomes the lowest energy state. 
The equivalent of this procedure in superconductivity is known as field cooling. 

Cooling under rotation has to be used in ■^He-A in order to stabilize and identify 
the single-quantum vortex (with n = 1) which at low flow velocities has lower energy 
than the double-quantum vortex which, in turn, has a much lower critical velocity. 
Again, the large difference in critical velocities of these two vortex structures arises 
because of their different core structures. As distinct from the doubly quantized vortex, 
the singly quantized ^He-A vortex has a hard vortex core with a radius comparable 
to the superfluid coherence length ^ (which lies hidden and embedded within a three 
orders of magnitude larger soft core of continuous structure |SH]). As a result, its 
critical velocity Vc ~ n/S, is close to the critical velocity for the formation of a B- 
phase vortex. The formation of the purely continuous texture of a doubly-quantized 
A-phase vortex does not involve breaking the superfluid state anywhere; it merely 
requires reorienting the degeneracy variables of the order parameter. That is why the 
corresponding critical velocity is much smaller, Vc ~ k/^d- 

These theoretical predictions have been tested in numerous different rotating 
experiments. Such measurements also indicate a wide range of variation in the 
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observed A-phase critical velocities, proving that they depend on the prehistory of 
sample preparation and thus on the quality of the global order parameter texture 
|59l IfiOj . Entry into the texture at unusually low critical velocity is provided by such 
regions where the spin-orbit coupling is not minimized, if they extend to the cylindrical 
boundary. Thus the largest reductions in critical velocity are observed in the presence 
of different types of planar domain-wall-like defects in the A-phase order-parameter 
field, which are called solitons If the plane of the soliton is oriented parallel to 
the rotation axis (called a splay soliton), then the critical velocity approaches zero, 
although it apparently always remains finite, and all vortex quanta will enter from 
along one of the two connection lines of the soliton sheet with the cylinder wall. This 
means that the emerging vortices are not lines but a periodic chain of circular and 
hyperbolic vortex quanta stacked within a soliton sheet (Fig. . This structure is 
called a vortex sheet It has the lowest critical velocity of all A-phase vortex 

structures and also the fastest dynamic response [021 ■ Numerical calculations of the 
flow instability of various different one-dimensional initial textures j63) show that semi- 
quantitative agreement with the measured variation exists. A further peculiarity of 
A-phase vortex textures is the existence of remnant vorticity in the form of vortex 
lines pinned in soliton sheets in the bulk liquid |()4| . 

The large difference between the critical velocities for the formation of a doubly 
quantized A-phase vortex and a singly quantized B-phase vortex makes it possible to 
create different metastable flow states in a rotating two-phase sample. 

2.4- Vortices and AB interface in rotation 

With two ■^He superfluids which belong to the same order parameter manifold it 
becomes possible to construct a unique situation which does not exist in other known 
coherent quantum systems. With a profiled magnetic-field distribution it is possible to 
achieve the coexistence of ^He-A and ■^Hc-B, and to stabilize an AB phase boundary 
in the sample. What happens to such a superfiuid two-phase sample in rotation? 
The substantial mismatch in the vortex properties of the two phases - their critical 
velocity, quantization of circulation, and vortex structure - raises the question: How 
are the vortices going to behave at the AB interface? In '^He-A the critical velocity 
is low, while in '^He-B it is an order of magnitude higher. This means that the A 
phase tends to be filled with essentially the equilibrium number of double-quantum 
vortex lines, while in the B phase there would at least initially be no vortices. Is such 
a situation stable, how is it going to evolve, and how are the single-quantum vortices 
of '^He-B going to fit in this picture if they emerge later at higher velocities? 

The left-hand side of Fig. 0] depicts the situation where the two-phase system 
is brought into rotation at constant temperature. When the rotation is started, A- 
phase double-quantum vortices are created at low critical velocity while no vortices 
are formed in the B phase. This expectation, confirmed by measurement |f)5| . means 
that the A-phase vorticity is not able to cross the AB interface and is accumulated on 
the A-phase side of the interface such that it coats the interface with a dense vortex 
layer. The layer is made up of a continuous texture of vorticity '66' and sustains the 
tangential discontinuity in the flow velocities of the superfiuid fractions on the different 
sides of the AB interface. Thus we have constructed a metastable state in which the 
two superfluids slide with respect to each other with a large shear-flow discontinuity, 
since the superfiuid fraction in the A phase rotates solid-body-like while that in the 
B phase remains stationary in the inertial frame. 
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Figure 4. Sketches of the AB interface in rotation. Left: In the A phase 
double-quantum vortex lines are formed at low rotation and we may assume 
it to be approximately in the equilibrium vortex state. At the AB interface 
the double-quantum vortices curve over to the cylinder wall and cover the AB 
interface as a vortex layer. This layer supports the discontinuity in the tangential 
velocities of the superfluid fractions at the AB interface. The width of the AB 
interface is on the order of the superfluid coherence length ^ while that of the 
vortex layer is three orders of magnitude larger, namely the dipolar healing length 
^D- Thus in this metastable state below the critical velocity the continuous A- 
phase vorticity does not penetrate through the AB interface, but gives rise to 
the unusual axial distribution in the flow velocities of the normal and superfluid 
components, i;n('') = S7r and Vs{r), as shown on the right for r = R. Right: In 
equilibrium rotation, an A-phase double-quantum vortex with winding number 
n = 2, dissociates at the AB interface into its 2tt constituents, the circular and 
hyperbolic merons, each with n = 1. Each meron gives rise to a singular point 
defect, a boojum, on the AB interface. The boojum is required as a termination 
point of a singular 2it (n = 1) B-phase vortex. Thus in the equilibrium state the 
continuous vorticity crosses the AB interface, transforming to singular vorticity. 
However, neither point or line singularities are easily created in superfluid -^116 
and therefore the vortex crossing takes place in a Kelvin-Helmholtz instability 
event. 



The minimum-energy state is shown on the right in Fig. ^ Here the vorticity is 
conserved on crossing the interface and both phases contain the equihbrium number of 
vortices at any given angular velocity of rotation. Accordingly, the number of double- 
quantum vortices in the A phase is one half of the number of singular singly-quantized 
B-phase vortices. On approaching the interface, the continuous A-phase vortex splits 
into its two 27r constituents or merons. Each of the merons ends in a boojum on the 
AB interface. The boojum is a point-like topological singularity of the orbital 1 vector 
at the AB interface, the termination point of a singular B-phase vortex on the AB 
interface. 

In practice singularities are not easily created in superfluid "^He: like in the case 
of vortices the energy barrier is too high compared to ambient temperature (typically 
by 6 - 9 orders of magnitude). Therefore the equilibrium state at the AB interface is 
not obtained by increasing rotation at constant temperature and pressure [SS] . Nor is 
it formed by cooling through Tc at constant rotation and pressure, which is the usual 
method to create the equilibrium vortex state below a second-order phase transition. 
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The reason is that here the first order A— +B transition also has to be traversed, to 
form the AB interface within the sample [HZIEEI- The closest approximation to the 
equilibrium state is obtained by starting with the equilibrium number of B-phase 
vortices in a single-phase sample at high rotation, with the barrier field at zero or at 
sufficiently low value. Next the barrier field is swept up (at constant fl, T, and P), 
until A phase and the AB interface is formed. In this case the equilibrium superfluid 
circulation is already trapped in the sample and cannot all escape. Finally, by reducing 
rotation to the point where B-phase vortices start to annihilate, one has reached the 
equilibrium vortex state. 

In contrast, if one simply starts increasing rotation of a two-phase sample with 
an AB interface, very different behavior is observed because of the energy barriers 
preventing the nucleation of point and line singularities. This is one of the ways 
to demonstrate the superfluid Kelvin-Hclmholtz instability, where the AB interface 
becomes unstable in the presence of an excessively large counterflow velocity which 
is oriented parallel to the AB interface (Sec. |3Jl. A complex chain of events is then 
started in which also vortices escape across the AB interface from the A to the B-phase 
side. 

2. 5. Vortex dynamics and mutual friction 

In •^He-B and ''He-II, in the absence of vorticity, the superflow is potential and the 
superflow velocity Vg = V"I> is expressed in terms of the "flow potential" which is 
proportional to the phase of the superfluid wave function, $ — {hl/2'k)lp. The superflow 
velocity Vg obeys the Euler equation |(i9j : 

^ + Vm = 0, (6) 

where jl — fi + and /i is the chemical potential. When quantized vortices 

(or continuous vorticity in "^He-A) are present, the superfluid velocity is no longer 
potential. The motion of a vortex leads to a phase-slip effect and modifles the r.h.s. 
of Eq. ©. 

Vortex lines form a part of the superfluid component, but the normal component 
influences their movement through mutual friction which arises from the scattering 
of normal excitations from the vortex cores. In the zero-temperature limit, where 
the normal excitations vanish, the motion of a vortex line is governed by the Magnus 
force only, so that the vortex velocity coincides with the local superfluid velocity at 
the position of the vortex element. At non-zero temperatures the friction between the 
vortex and the normal component - the so-called mutual friction - causes a drag force 
on the vortex line and, as a result, the velocity vl = dvi^/dt of a vortex segment in 
the flow deviates from Vg. 

In the presence of vortices, the flow potential is not uniquely deflned along the 
contours encircling the singular vortex lines. If the vortices do not overlap, the flow 
potential can be written as 

N 
/3=1 

Here Tp{sp) are the coordinates of singular lines specifled by a parameter sp. If the 
positions of these lines also depend on time, = rp{t), the time derivative of the 
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superflow velocity becomes 



/3 ^ 



Here is the derivative only of the explicit t dependence of $ and we put 

V/3 — dvp/dt. The vorticity of a single vortex is 

ojf) ~ curlvg^ = curl V<i>^ . 

In ^He-B and ''He-II the vorticity from singular vortex lines is expressed as 

u}j3^ Kp J ^i^- ^p) drp . (7) 

Here rp is the coordinate of the /3-th vortex line and S {r — rp) is the three-dimensional 
(5-function, drp = sp dsp, sp is the unit vector in the direction of the vortex hue at the 
point Tp, and dsp is the arc length of the vortex line. The circulation of each vortex 
Kp = npK may have np circulation quanta k. Since the derivative 

p 

is usually defined as the "superfluid chemical potential" /ig = —d^/dt we have 
'dt 



^ + Vms = ^v^ xw^. (8) 

P 

One can write here fis = A + Ms, where jis — Ms — M is the deviation of the 
superfluid chemical potential from the total chemical potential of the fluid, which is 
the counterpart of the gauge-invariant scalar potential in the theory of nonstationary 
superconductivity j7()| . In equilibrium fis = 0, so that /is — M, which in fact is the 
famous Josephson relation. 

The velocity of each vortex is determined up to its component perpendicular to 
the vortex line |27|: 

vp = sp X (vs X S;3) + asp X (v„ - Vg) - a'sp x [sp x (v„ - v^)] (9) 

Here a{T,P) > and a'{T,P) < 1 are the temperature and pressure dependent 
dissipative and reactive mutual-friction parameters. 
Inserting Eq. © into Eq. JHJ we find 

^vp X LOp ^ ^Kp S{r~ rp) dsp [sp x [sp x {sp x v^)]] 

P P 

+ a'Y^Kp I 5{y - rp) dsp [sp x [sp x {sp x (vn - v^))]] 

P 

- a^^Kp 5{r-rp) dsp [sp x [sp x (vn - v^)]] . 

P 
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The first line gives 

6(r - dsp [sp X [s^ X (s^ x v^)]] 

13 

^ ^ X! '^'3 M(r - r/3) dsf) [sf) X Vs] = Vs x a;s 

13 

where 

a;s = ^a;/3 (10) 

is the total vorticity of the superfluid. The second line can be transformed in the same 
way. As a result, 

+ V/i = Vs X a;s + fmf . (11) 
Here fmf is the mutual-friction force 

fmf = - a ^ / (5 (r - r/3) C?I",3 X [s/5 X (Vn - Vg)] 

+ a' [(vn - V,) X a;,] (12) 

exerted by the normal component on a unit mass of superfluid via the vortex lines. 
The first term in fmf is the viscous component, with a negative projection on the 
relative velocity Vs — Vn, 

fmf -(vs-Vn) = -(vs-Vn)^a^K^ M (i" " r^) [1 - cos^ 7^] rfs^.(13) 

13 •' 

Here 7/3 is the angle between and Vg — Vn. Without mutual friction force, Eq. (fTT|l 
coincides with the Euler equation of the classical hydrodynamics of an ideal inviscid 
fluid 

d'v 

— + V/i = vxu;. (14) 
at 

In what follows we consider the case where all vortices have the same circulation with 
n=l. 

For '^He-B, the mutual-friction parameters were measured in Refs. |721IZS|- This 
remarkable hydrodynamic measurement was performed by examining capacitively 
changes in antisymmetrically driven normal modes of a thin membrane which lies 
in the plane perpendicular to the rotation axis. As will be seen in Sec. [S] (Fig. I2ti|) . 
these results agree with later measurements on the longitudinal propagation velocity 
of vortices in the rotating cylinder which also yield the dissipative parameter a. These 
experimental data agree with calculations of the parameter values in Refs. (for 
a theory review see |11|'). The mutual friction parameters are discussed also in terms of 
the chiral anomaly and the Callan-Harvey effect in relativistic quantum field theory in 
Ref. 0]. Both a and a' vanish at T = 0. With increasing temperature the dissipative 
mutual-friction parameter a increases so that above 0.6 all dynamic processes are 
heavily overdamped in ^He-B. In '^He-A a is expected to be in the overdamped regime 
at all currently accessible temperatures P I72L 171] . In contrast, in ''He-II a is much 
smaller at all experimentally relevant temperatures, as seen in Fig. |31 The difference 
in the magnitudes of the mutual- friction parameters for superfiuid ^He and "^He is 
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T/T^. T(K) 

Figure 5. Mutual friction covers a different range of values in '^He-B (left 
column of panels) and ''He-II (right column). The top row of panels shows 
the dissipative mutual friction coefficient a{T) and the bottom row the dynamic 
parameter = (1 — a')/a. A transition in vortex dynamics is expected at 
g ~ 1. The data for ^He-B is from Refs. EUESl and for "^He-II from Ref. 



SO striking that the usual picture of '^He-II vortices moving with the superflow fails 
for '^He vortices which move across the flow in most of the accessible temperature 
range. In this respect, motion of vortices in '^He resembles that in most of type-II 
superconductors where they also move perpendicular to the flow with practically no 
measurable Hall angle [Ilj . 

Since the equation of motion for the trajectory of a vortex only depends on 
the two parameters a and a', its solutions can be classified in terms of the dynamic 
parameter q defined as 



This parameter has already been introduced in Section 11.11 According to both the 
theoretical predictions jJOl and the experimental data for ■^He-B the parameter q 
depends exponentially on temperature in the low-temperature limit. The lower row of 
panels in Fig.[Slshows the inverse = {l — a')/a for ■^He-B and "'He-II. In superfluid 
dynamics q acquires an important function: it characterizes the relative influence of 
dissipation and 1/q plays the role of the Reynolds number fSec. 14.31) . ^ change in 
the characteristic solutions and a corresponding transition in the dynamics can be 
expected in the regime q ^ 1: this is located in the middle of the temperature range 
for '^He-B, but only a few /iK below T\ for ''Hc-II. 

It is instructive to inspect some solutions of Eq. as sketched in Fig. El for an 
ideal rotating cylinder where the rotation and cylinder axes coincide. We assume the 
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experimentally important situation where the applied counterflow velocity dominates, 
i.e. we set in Eq. jnj Vn = O x r and Vg « 0, and consider a large enough q so that the 
vortex dynamics is regular (instead of turbulent, see Sec.^. The simplest case (a) is 
a rectilinear vortex line parallel to the rotation axis which is released in the flow from 
the cylinder wall. This might correspond to the situation in which a double-quantum 
vortex forms in ^He-A when the cylinder is filled with a perfect global texture. In the 
cylindrically symmetric flow geometry the vortex will retain its rectilinear shape while 
it travels inward along a spiral which finally places the line in its equilibrium position 
along the axis of the cylinder. Another important case is (b) a vortex with a free end 
expanding in a long cylinder. The free end terminates perpendicular on the cylinder 
wall and describes a regular helix on the wall during its expanding motion. If we 
neglect the influence from its own curvature on the motion, then the curved section of 
the vortex remains within the same radial plane which rotates at the angular velocity 
(1 — a')fl around the cylinder axis, when viewed from the rotating frame. The trailing 
section of the vortex becomes aligned along the central axis and remains at rest there. 
An immediate extension of this example is case (d) where the vortex expands in a 
cylinder in the presence of a pre-existing cental vortex cluster. The rectilinear trailing 
end of the vortex becomes now part of the vortex cluster and is incorporated as one of 
the vortices in the peripheral circle of lines. This section would prefer to be stationary 
in the rotating frame, it resides in a region of the sample cross section where the 
average counterflow velocity vanishes, and therefore it moves only with difficulty from 
one lattice site to the next in the outermost ring of vortex lines. The free end, on the 
other hand, expands along a spiral trajectory and leaves behind a helix which cannot 
relax instantaneously. 

To explore the hydrodynamic transition between regular and turbulent vortex 
dynamics at g ^ 1, we can ask the following question: how is the dynamic equation 
Hll|l modified for superfluids in the continuous limit, after averaging locally over vortex 
lines? This can be done even for tangled vortex states, if the lines are locally sufficiently 
parallel and their radius of curvature is much larger than the vortex core diameter. 
This would be the case, for instance, for rotating states with transient time-dependent 
disorder, since the rotating flow would rapidly polarize the component parallel to the 
rotation axis while the order in the transverse plane would be restored only later. In 
this case the averaging over nearby vortices on the r.h.s. of Eq. gives 

dv 

— + V/i = Vs X u;s-a'(vs-Vn) x uj^+a lj^ x [Ws x (v^ - Vn)] .(16) 
ot 

This result is known as the coarse-grained hydrodynamic equation for the superflow 
velocity [75]. It will later be used to develop an analogous superfluid interpretation 
to the transition to turbulence, as can be derived in viscous hydrodynamics from the 
Navier-Stokes equation |25|. 

2.6. Kelvin-wave instability of vortex lines 

A quantized vortex can support helical Kelvin waves which become important in 
turbulent vortex dynamics. The Kelvin waves may lead to an instability which 
generally develops at large flows. This instability can result in an increase in the 
number of individual vortex lines under an applied flow if the mutual friction is 
sufficiently small. It is the first step in the process by which more vortices can be 
generated from one or a few existing vortices so that bulk volume turbulence can switch 
on. The evolution towards the instability starts when a vortex becomes sufficiently 
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a) rectilinear b) curved end c) lialf ring d) curved end outside 

vortex of single vortex on cylinder wall vortex cluster 



Figure 6. Sketches of vortex trajectories in an ideal rotating cylindrical superfluid 
sample, as viewed in the rotating frame, (a) A rectilinear vortex, which is released 
in the flow from the cylinder wall, travels along a spiral to its equilibrium location, 
to become aligned along the center axis of the cylinder, (b) In a long cylinder 
the propagating end of the vortex describes a spiral trajectory on the cylinder 
wall while it expands into the vortex-free section of the cylinder. The trailing 
end is aligned along the central axis, the equilibrium location for the vortex. 
(c) An elementary vortex usually forms as a half ring on the cylinder wall. It 
expands both axially and radially. Neglecting the contribution from its self- 
induced velocity, the half ring positions itself perpendicular to the azimuthally 
circulating applied flow and remains during its expansion inside a radial plane 
which rotates around the cylinder axis with angular velocity (1 — Q!')r2 with respect 
to the cylinder walls, (d) If a central vortex cluster already exists in the cylinder 
(with A'^ = ttR^ 2Q/k rectilinear vortex lines) , then the spiralling motion leaves 
behind a helically wound vortex. At low temperatures, when the axial motion 
is slow, a tightly wound helix forms around the cluster. The end section of this 
vortex, where it approaches the cylindrical wall, is most likely to suffer the Kelvin- 
wave instability. Here the applied flow velocity is large and the vortex is partly 
aligned azimuthally along the applied flow. 



aligned along the applied flow. A vortex line oriented parallel to the external flow is 
an unstable configuration, as demonstrated by Donnelly et al. by applying a thermal 
counterflow current parallel to rectilinear vortex lines in rotating '*He-II 1 77J . Above a 
critical axial flow velocity rectilinear vortices became unstable and formed a turbulent 
tangle with varying axial polarization, depending on the axial counterflow velocity. 
The phenomenon was explained by Glaberson et al. |78l I79| in terms of the Kelvin- 
wave instability and was later reproduced also numerically by Tsubota et al. |8()| . 

Consider small transverse deformations to an isolated vortex line, oriented along 
the z axis, in externally applied counterflow v — v z parallel to it. Parametrizing the 
position vector of an arbitrary element on the deformed vortex as 

r{z,t) = Ciz,t)x + Tjiz,t)y + zz, (17) 

we can write the unit tangent of the line as s(z) « {d(/dz)x + {dr]/dz)y + z (to 
linear order in the small quantities rj). The vortex curvature also gives rise to a 
self-induced contribution to the superfluid velocity at the vortex line. In the local 
approximation [23 ESI 

where k = (K/47r) ln(27r/fcrc). Inserting these to the equation of motion for the vortex 
Hne, Eq. ©, we flnd 

dC, d'^rj ( d'^C dri\ , ( d'^rj dC, 

— = - kttt + oi K-—1T + V— + a K—^ - V— 
at oz'^ V oz'^ oz I \ az^ az 
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) 



(19) 



The dispersion relation for wavelike disturbances oc exp[— ^(tjt — kz)], i.e. Kelvin 
waves, can be found as 



Here we again encounter the dynamic parameter q — a/(l — a') already introduced in 
Eq. (fT5|l as the ratio of the imaginary and real parts of the dispersion |H2| . The waves 
are always damped but at high temperatures, where q > 1, the waves are overdamped 
and do not propagate. 

On the other hand, in the presence of externally applied flow the long-wavelength 
modes with k < v/h have Im(a;) > 0, and exhibit exponential growth. In other words, 
if an evolving vortex configuration at some time has long enough vortex-line sections 
oriented parallel to external flow, these will become unstable to exponentially growing 
helical deformations. The expanding waves can then undergo reconnections, either 
with the walls of the container or with other vortex lines. This leads to a growing 
number and density of vortices and, ultimately, if q is small enough, to the onset of 
turbulence, as discussed in Sec. 14.61 

3. Kelvin-Helmholtz instability in superfluids 

3.1. Introduction 

The Kelvin-Helmholtz shear flow instability is a well-known phenomenon of classical 
hydrodynamics which was first discussed by Lord Kelvin already in the 1860's. It 
occurs at the interface between two fluid layers which are in relative motion with 
respect to each other. For instance, at low differential flow velocity the interface 
between two stratified layers of different salinity or temperature is smooth in the 
ocean, but at some critical velocity waves are formed on the interface. Similarly, 
ripples do not form on the water surface on a lake at infinitesimal wind velocities, but 
form at some finite critical value. The Kelvin-Helmholtz instability is thus a common 
phenomenon in nature around us. The condition for instability is derived in many 
textbooks on hydrodynamics j8HI I84j in the limit of inviscid and incompressible fluids. 
In superfluids it was first observed at the rotating AB interface |65| . 

The AB interface instability occurs between the two ^He superfluids, '^He-A and 
^He-B, when their superfluid fractions move with respect to the normal component 
tangential to the AB interface (Fig.0]/e/t). This initial state is dissipationless while the 
state after the instability is not, as surface waves or ripplons form on the interface and 
their motion is highly damped. In conventional liquids and gases the mathematical 
description of the interfacial instability is inevitably only approximate: if viscosity is 
neglected, the initial nondissipative states of the two liquid or gas layers sliding with 
respect to each other are not exact. The question then arises not only about the true 
value of the critical velocity, but more generally about the existence and nature of the 
instability. Superfluids are the only laboratory examples of cases where viscosity is 
totally absent, and the mathematical description of the instability can be presented 
analytically in a simple form. The initial state of the AB interface with different 



cj = (1 — a')kk'^ — a'kv — ia{kk'^ — kv). 
With vanishing counterflow, w = 0, the dispersion relation simplifies to [5T] 
Lo — — a' — ia). 



(20) 



(21) 



Dynamics of vortices and interfaces in superBuid ^He 



24 



tangential superflow velocities across the interface is in nondissipative thermodynamic 
equilibrium, until the critical velocity of the hydrodynamic instability is reached. As 
a result - contrary to viscous normal fluids - the instability threshold is well defined. 
Experimentally it is manifested by the sudden formation of vortices in the initially 
vortex- free ^He-B phase (Fig.^right) at a rotation velocity at which no vortices would 
yet be formed without the presence of the AB interface. 



3.2. Kelvin- Helmholtz instability in classical hydrodynamics 

The Kelvin-Helmholtz (KH) instability is one of the many interfacial instabilities in 
the hydrodynamics of liquids, gases, charged plasma, and even granular materials. 
It refers to the dynamic instability of an interface with discontinuous tangential flow 
velocities and can loosely be defined as the instability of a vortex sheet. Many natural 
phenomena have been attributed to this instability. The most familiar ones are the 
generation of capillary waves on the surface of water, first analyzed by Lord Kelvin 
|85|. and the flapping of sails and flags, flrst discussed by Lord Rayleigh |86| . 

Many of the leading ideas in the theory of interfacial instabilities in hydrody- 
namics were originally inspired by considerations about ideal inviscid flow. The 
corrugation instability of the interface between two immiscible ideal liquids, with 
different mass densities pi and p2, occurs in the gravitational field at the critical 
differential flow velocity 

(Ui - U2) = 4(T.g , (22) 

P1P2 

where a is the surface tension of the interface, and g is the gravitational acceleration. 
To separate the gravitational and inertial properties of the liquids, we rewrite the 
threshold velocity in the form 
P1P2 



-iVi-V-,y = 2VaF . (23) 

pi + P2 

Here F is the external held stabilizing the position of the interface, which in the case 
of a gravitational field is 

F = gipi - P2) , (24) 

but can in general originate from some other source. The surface mode of ripplons, 
or capillary waves, which is first excited has the wave number corresponding to the 
inverse 'capillary length' 



fco = VF/^ ■ (25) 

However, ordinary fluids are not ideal and the correspondence between this theory 
and experiment is not satisfactory. One reason for this is that the initial state cannot 
be properly prepared - the shear-flow discontinuity is never an equilibrium state in 
a viscous fluid, since it is not a solution of the Navier-Stokes equation. As usual in 
hydrodynamics, it is not apparent whether the notion of 'instability' can be properly 
extended from the idealized inviscid model to the proper zero viscosity limit of the 
Navier-Stokes equation and further to the case of finite viscosity [S7| . 
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3.3. Experimental setup 

The superfluid Kelvin-Helmholtz experiments reported in Ref. i6.5J were conducted in 
a rotating nuclear demagnetization cryostat in which the hquid '^He sample can be 
cooled below 0.2 Tc in rotation up to 4 rad/s. A schematic illustration of the sample 
setup is shown in Fig. 13 The sample is contained in a fused quartz tube of i? = 3 mm 
radius. An aperture of 0.5 - 0.75 mm diameter in the bottom end plate restricts the 
flow of vortices into the sample from the heat exchanger volume below. The quality of 
the sample-tube surface is of great importance, since any larger surface defects serve 
as nucleation and pinning sites. Before use the container is treated with a HF solution, 
to etch away sharp protuberances, and then cleaned with solvents. The same quartz 
tube may display sizable variation in critical velocity and vortex formation properties 
depending on how it is cleaned between experiments, indicating that dust particles, 
dirt, frozen gas, etc. serve as nucleation and pinning sites. Most measurements were 
performed using a sample tube which allowed vortcx-free rotation up to the maximum 
rotation velocity of the cryostat. 

Two independent continuous-wave NMR spectrometers are used for monitoring 
the sample. Their detector coils are located close to both ends of the sample volume. 
In these regions of the sample a low magnetic field is needed for NMR polarization 
which has to be as homogeneous as possible. Therefore three separate superconducting 
solenoidal magnets with compensation sections, all with their fields oriented along 
the rotation axis, are needed in the experiment. Two magnets are required for the 
homogeneous NMR fields and one as a "barrier magnet" around the central section 
of the tube for stabilizing ^He-A. The barrier magnet provides fields up to 0.6 T 
which is sufficient to stabilize the A phase at all temperatures and pressures. In the 
inhomogeneous barrier field H{z) the equilibrium position of the AB interface is 
determined by the equation H{zq) — iJAB(r), where Hab{T) is the first-order phase 
transition line between the A and B phases in the {H, T) plane (Fig. |H1 left) . 

The two superfluids, ^He-A and ^He-B, have very different magnetic properties, 
as seen in the phase diagram on the left of Fig. |S1 In the region where H > Hab (T) 
the A-phase has lower magnetic energy, while in the neighboring region H < Hab (T) 
the B-phase is favored. Thus the gradient in the magnetic energy densities of the two 
liquids provides a restoring force F on the AB interface, 



Here XA > Xb sue the magnetic susceptibilities of the A and B phases, respectively. 
Accordingly, this VH dependent restoring force on the AB interface has to be used 
instead of gravity in the KH instability condition. 

Various configurations of A and B phases can be stabilized and trapped in the 
sample tube, (i) At zero or low barrier fields the B phase might occupy the whole 
sample volume. In this all-B configuration the two spectrometers probe the same 
volume. The time difference between the two NMR readings can then be used to 
study the propagation of vortices along the column, (ii) At barrier fields above 
Hab {T, P) , the equilibrium transition field between the A and B phases, the A phase 
is stabilized in the center of the barrier magnet. In this BAB configuration the top and 
bottom sections of the sample are disconnected: vortex lines generated in one of the 
B-phase sections do not pass to the other across the A phase region, (iii) At pressures 
above 21 bar, ^He-A is stable also at zero magnetic field from Tc down to Tab{P)- 
At high temperatures the entire sample volume is then filled with A phase (all- A 




2 V [(XA - Xb)H'] 



(26) 
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Figure 7. Two-phase superfluid ^He sample for measurements on the Kelvin- 
Helmholtz instability. The sample volume is 6 mm in diameter and 11 cm long. 
A small orifice of 0.5 - 0.75 mm in diameter in the center of the bottom end plate 
separates the sample from the rest of the liquid ^He volume. The ''He below the 
orifice is needed to establish thermal contact to the refrigerator. NMR pick-up 
coils are located at both ends of the sample tube. They are circular coils with 
their symmetry axis transverse to the tube axis. Two solenoidal superconducting 
magnets provide the homogeneous axially oriented polarizing fields for NMR. A 
third barrier magnet creates the field to stabilize ^He-A in the center section of the 
long sample. Two examples of the A phase region (shaded) are shown at 0.55 Tc'. 
at the current of /b = 8 A in the barrier solenoid the A phase extends further in 
the column and the two AB boundaries are almost flat, while at = A A the A 
phase region is short and the AB boundaries are concave. 



configuration). When the temperature is lowered somewhat below the equilibrium 
A— s-B transition, B phase nucleates in the heat-exchanger volume and expands into 
the lower section of the sample (at about 0.75 Tc)- If the barrier field is sufficiently 
high to stabilize the A phase, the advancement of the B phase is stopped by the stable 
A phase and an AB interface is formed. Since the A phase in the top section can 
supercool quite substantially '89J , the sample remains in the BA configuration (Fig. |S1 
right). Thus the barrier field isolates the top from an A— >B transition in the lower 
section. Eventually, at a low enough temperature, B phase also forms independently at 
the top. Since the ultimate supercooling of "^He-A depends on the surface properties, 
the A— >B transition also serves as a measure of the quality of the quartz walls. At 
best the top was supercooled to 0.52 29 bar (when the equilibrium transition 

is at Tab = 0.85 Tc)- The BA configuration was most important in the early stages 
of measurements on the AB interface instability, since the top spectrometer is then 
recording the vortex number in the A phase and the bottom spectrometer in the B 
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Figure 8. Magnetic stabilization of '^He-A in the experimental setup of Fig. 171 
1651 . (Right) Sample container and the axially oriented magnetic field distribution 
along the sample in the BA configuration (with current 7^, = 4.0 A in barrier 
magnet). The barrier field maintains A phase in stable state in the section in 
the middle of the sample, where H > Has{T). In the top and bottom sections 
only the NMR polarization field is applied and H < H\q{T). In the top A phase 
may persist in meta-stable state (see text for explanation). One AB interface 
exists here at the lower location H{—\zo\) = HxsiT), where z = —\zq\ is the 
average axial location of the interface. When the temperature T or the current 
7b is varied, changes occur in both zo(T, /b) and in the field gradient Vi7 in 
which the interface resides. The restoring force acting on the AB interface in 
Eq. 1261 depends on the field gradient Vi7 and therefore on the field profile of 
the barrier magnet. (Left) Magnetic phase diagram in the (H,T) plane, where 
Hab{T) marks the first-order phase transition line between A and B phases (at 
P = 29.0 bar pressure) 1881 . The upper plane at 20(^1 ^b) separates the section in 
the middle, where '^He-A is the true vacuum state, from the section in the top, 
where ^He-A is the false vacuum. The false vacuum persists down to a container 
dependent minimum temperature. 



phase. To understand what happens in the two-phase sample when rotation is started, 
and to correlate in real time the events on both sides of the AB interface, simultaneous 
recordings of vortex numbers in both phases are needed. 

It turned out that in rotation '^He-A is in solid-body-like rotational flow, 
practically locked to co-rotation with the container, owing to its low critical velocity of 
vortex formation. In contrast, in '^He-B the superfluid component is in the vortex-free 
Landau state and stationary in the laboratory frame. This nondissipative initial state 
becomes possible through the formation of an A-phase vortex layer which covers the 
AB interface and provides the discontinuity in the tangential flow velocities (Fig. 01 
left) . The tangential discontinuity is ideal - there is no viscosity in the motion of the 
two superfluids so that this state can persist for ever. If the rotation velocity is then 
incrementally increased to the critical value, the Kelvin-Helmholtz instability occurs 
and some vortices penetrate across the AB interface into the vortex-free B phase (Fig.^] 
right). The first time, when this happens in increasing rotation, is called the critical 
rotation velocity f2c(T', ^b)- If rotation is increased further, then the instability will 
occur repeatedly every time when the counterflow velocity close to the cylinder wall 
reaches the critical value. In each case, the signal for the instability is the appearance 
of new vortices in the B-phase section. 

One might think that by substituting the interfacial restoring force F ()26|l into 
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Eq. H23|l . and using the superfluid densities of the A and B phases instead of the 
total density, the critical velocity for the KH instability of the AB interface could 
be obtained. However, it turns out that a proper extension of the KH instability to 
superfluids incorporates the criterion in Eq. H23|l only as a particular limiting case, 
namely when T — > 0, while in general the physics of the instability is rather different 
from the ideal inviscid model. A different but well-determined criterion is obtained 
for the instability condition in terms of the velocities and densities of the superfluid 
fractions. This criterion is also formulated in the absence of viscosity, as one would 
expect in a perfect superfluid environment. Therefore the tangential discontinuity at 
the interface between '^He-A and "^He-B below the instability is a stable state. 

3.4- Modification of Kelvin- Helmholtz instability in superfluids 

The criterion for the KH instability of ideal fluids in Eq. Il23(l depends only on the 
relative velocity across the interface. In practice there always exists a preferred 
reference frame imposed by the environment. In the superfluid case this is the frame 
fixed to the container. At T ^ the normal component provides the coupling to the 
reference frame: in the state of thermodynamic equilibrium the normal component 
moves together with the container walls, i.e. Vn = in the frame of the container. 
Owing to this interaction of the AB interface with its environment the measured 
instability occurs at a lower differential flow velocity, before the classical criterion in 
Eq. H23() is reached. Moreover, it even occurs in the case when the two superfluids 
move with the same velocity. However, this does not imply that the renormalized 
instability criterion would depend on the interaction with the normal component - 
in fact, it is still determined only by thermodynamics. Waves are formed on the 
interface when the free energy of a corrugation becomes negative in the reference 
frame of the environment and the AB interface becomes thermodynamically unstable, 
i. e. when it becomes possible to reduce the energy via the normal component and its 
interaction with walls. The instability condition Eq. (|23|l of the ideal inviscid fluid is 
restored if the interaction with the environment is not effective. This would occur, 
for example, during rapid rotational acceleration at very low temperatures when the 
instability caused by the interaction with the environment does not have enough time 
to develop. 

The free energy of the disturbed AB interface C{x, y) in the reference frame of 
the container contains the surface tension energy, which corresponds to the potential 
energy in the 'gravity' field, and the kinetic energy of the two liquids 



Here we take into account that in thermal equilibrium the normal component is at 
rest in the container frame, — VnB = 0, and only the superflow contributes to the 
kinetic energy. The heights of the A and B phase layers are denoted with /ia and h-Q. 
For simplicity we ignore the anisotropy both in A phase and in B phase at high fields 
and approximate the superfluid-density tensor with a scalar. We write the superfiuid 
velocities as VsB(r) = Ub + <5vsb(i") and VsA(r) = Ua + (5vsA(r), where Ub and Ua 
are the velocities at an undisturbed flat interface. Using V x Vg = V • Vg = 0, and the 
boundary conditions at the interface, s • VgA = s • Vsb = 0, one obtains the free energy 




(27) 
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Figure 9. AB interface instability as a function of temperature. The measured 
critical velocity flc = Ub/R is plotted at three different fixed currents I^, in the 
barrier magnet Ib5l . When T is changed at constant I^, the AB interface relocates 
roughly along the contour where H{z, r) = Hp^b{T). Thus its average position zq 
and the field gradient V-ff change. These changes alter the value of the restoring 
force in Eq. 1261 . The solid curves represent the instability criterion in Eq. l3Ui . 
when Ub = and Ua = ''nA = ''nB = 0. No fitting parameters are used. The 
vertical arrows below the figure demarcate the temperature interval where the AB 
interface exists within the barrier magnet at different currents At TABC^^min) 
the interface is stable in the minimum field location. At Tab (^^(^b)max) the 
critical field i^AB equals the maximum field inside the barrier magnet. 



of a surface mode as (for details see Ref. where also the anisotropy of the tensor 
Pg-J^ is taken into account) 

J^{Ck)^a^F + ea-k{pful + pfUi)] . (28) 

Here k is the wave vector along x for a surface corrugation amplitude of the form 
Cix) = asinkx, and 

cff ^ PiiA off ^ PsB .2q^ 

'^^ tanh(/c/iA) ' ^ tanh(fc/iB) ' 
In the relevant experimental conditions we are always in the 'deep water' limit, 
khA ^ 1, khB ^ 1- The free energy becomes negative for the first time for the 
critical ripplon wave number ko = (F/aY^^ when 

^ (PsbC/| + PsaUI) = . (30) 



r2 I „ tt2^ 

A comparison of Eq. (|30|) to the measured critical rotation velocity ilc of the first KH 
instability event is shown in Fig. |5| The curves have not been fitted; they have been 
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drawn using accepted values from the literature for the various quantities (for details 
see Ref. ^^)- Such a remarkable agreement for a complicated phenomenon can only 
be achieved in supcrfluid ^He! 

We thus conclude that, even under perfectly inviscid conditions, in superfluids 
the critical velocity for the KH instability is not given by the classical result for ideal 
fluids. The new criterion for the instability in Eq. (j^OI) has at first glance paradoxical 
consequences. The instability is not determined by the relative velocity Vsa — Vsb; in 
fact, the instability would occur even if the two superfluid velocities were equal. The 
instability would also occur for only a single superfluid with a free surface. These new 
features arise from the two-fluid nature of the supcrfluid. Therefore, the instability 
threshold is determined by the velocities Vg^ — Vn of each superfluid i with respect 
to the reference frame of the walls and thus with respect to the normal fractions 
which in thermodynamic equilibrium move together with the walls. Accordingly, the 
free surface of a superfluid - the interface between the superfluid and the vacuum 
- becomes unstable when the superfluid velocity exceeds the critical value in the 
reference frame of the normal fraction |22[ I23j . With many (i) superfluid fractions in 
the same liquid, such as the neutron and proton superfluids in a neutron star [021 ESI j 
the threshold is determined by some combination of the superfluid velocities Vgi — Vn 

m 



3.5. Kelvin- Helmholtz instability in the low-temperature limit 

On approaching the zero-temperature limit the density of normal excitations rarefles, 
the coupling with the container walls becomes weaker, and the superfluid density 
becomes the total density: psA PA and psB Pb- How is the superfluid going 
to react to the environment under these conditions? Let us compare the result in 
Eq. I|30|l with the ideal classical condition in Eq. (|23|l in the limit when T — 0. The 
classical instability condition reads 



(Ub - Ua)' - 2VaF , (31) 

PA + PB 

while the superfluid instability occurs when 

i {paUI + pbUI) = . (32) 

In the experiment the A-phase superfluid component is approximately stationary in 
the container frame, Ua ~ 0, and the densities of the two liquids are practically equal. 
Pa ~ Pb — P- Then the B-phase critical velocity at the instability is J/g = 4\/ uF j p 
according to the classical criterion, while from Eq. (|3()|l we obtain a result which is 
smaller by a factor of V2, i-e. Vb ~ 2\/ aF j p. 

The difference is imposed by the interaction with the environment 4 . To 
understand this, we repeat the derivation of the classical KH instability 94 with 
one important modiflcation. We add to the equation of motion of the AB interface a 
friction force which arises when the interface moves with respect to the environment, 
i.e. with respect to the normal component. In the reference frame of the container it 
has the form 

^ftiction = — r^tC • (33) 

In the low-T limit, the friction between the AB interface and the normal component 
arises from the Andreev scattering of ballistic quasiparticles from the interface. In 
this regime the parameter T oc [SSj EH E3 • For simplicity, we choose a situation 
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Figure 10. Sketch of the imaginary and real parts of the frequency uj{k) of the 
critical ripplon (when k = kg) at the interface between "^He-A and '^He-B. The 
diagram is constructed in the rotating frame with v^j\ = = and 

Pa ~ Pb = P on approaching the limit T — ► 0, and considering only incremental 
changes from the critical conditons. At the superfluid instability, the imaginary 
part Im uj{k) crosses zero as a function of tigB and the attenuation of ripplons 
transforms to amplification. Simultaneously the real part Re Lo{k) also crosses 
zero. The region where Re oj < and where the ripplon has negative energy is 
called ergoregion. The slope of the imaginary part is proportional to the friction 
parameter F. If T is strictly zero, and thus the connection with the environment 
is lost, the surface instability occurs at the value of v^b larger by a factor \/2 - 
at the corner point which is part of the branch obeying the ideal classical KH 
criterion of Eq. 1231 . 

when both fluids move along the x direction, and consider surface waves (ripplons) 
of the form C(^i^) = asui{kx — cut) in the container frame. The classical dispersion 
relation for ripplon motion is then modified by the presence of the friction term to 

/cj \2 /u \2 F + k'^cr uj , , 

Pa(^-C/a) +Pb(^-C/b) =^ .r-, (34) 

modifying the nature of the instability. We rewrite the above equation as follows: 
k PA + Pb 

± , , W 7 *r- {Uk-UbY . (35) 

VPA + Pb V k k PA + Pb 

If r = 0, this reduces to the classical dispersion relation and the instability occurs when 
the expression under the square root becomes negative. The ripplon spectrum then 
acquires an imaginary part with both plus and minus signs. Thus at the threshold 
velocity of Eq. H31(l the perturbation grows exponentially in time. A sketch of the 
imaginary and real parts of the frequency of the critical ripplon mode (with k = ko) 
is shown in Fig. IIUI 

In the case of the superfluid instability we have to consider the influence of the 
term with the friction parameter T in Eq. 1351 When the imaginary part Im Lo{k) 
crosses zero as a function of Ub, the attenuation of the ripplon modes is transformed 
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to amplification. This occurs first for ripplons with the wave vector given by Eq. H25|l . 
While the instability condition Eq. (|32|l does not depend on the friction parameter 
r, the slope of the imaginary part is proportional to T. Therefore, when T is strictly 
zero and the connection with the reference frame vanishes, the interface instability 
starts to develop when the classical KH criterion in Eq. H31|l is reached. In the 
limit T and F — > 0, the experimental result is expected to depend on how 
the observation time compares to the time needed for the interface to become coupled 
to the environment, and for the superfluid instability to develop at the lower critical 
velocity. With sufficiently rapid rotational acceleration the classical KH instability 
might then be reached. Unfortunately at present no experimental verification exists 
on these predictions. 

The real part of the ripplon frequency also crosses zero at the same velocity than 
the imaginary part. Above the threshold the real part of the ripplon spectrum, and 
thus its energy in the container frame, becomes negative. This creates a connection 
to the physics of black holes |HH1- In general relativity the region where a particle has 
negative energy is called the ergoregion. In the 'shallow water' limit kh <C 1, when 
the spectrum of ripplons becomes 'relativistic', an exact analogy with the relativistic 
quantum field in the presence of the black- or white-hole horizons emerges ^01 . 

It is also possible here to identify a similarity with the Kelvin-wave instability of 
an isolated vortex line in applied flow. In the T ~^ limit, when both a ^ and 
a' —^ 0, the instability of a vortex line in externally applied superflow towards the 
formation of Kelvin waves resembles the AB interface instability. For the unstable 
modes in Fig. EH with Im a; > 0, the real part of the spectrum is negative. Re w < 0. 
Similarly, for a vortex of finite length L the wave-number is constrained from below, 
k > ko = 2Tr/L, and the instability forms at a critical velocity Vc kk^. This Vc does 
not depend on the friction parameter a, whose role is to provide the dissipation from 
the interaction between the vortex and the environment (i.e. the normal component). 

To conclude, we have found that the central property of the superfluid KH 
instability is that the instability condition does not depend on the relative velocity 
of the superfluids, but on the velocity of each of the superfiuids with respect to 
the environment. The instability occurs even if the two fluids have equal densities, 
pA = Pb, and move with the same velocity, Ua = Ub- This situation resembles that of 
a flag flapping in wind. It was originally discussed with the view of the KH instability 
of ideal fluids by Lord Rayleigh [HS| ■ One might assume the instability to be that of a 
passive deformable membrane between two distinct parallel streams having the same 
density and velocity, as in the superfluid KH example: the flag being represented by 
the AB interface, and the flagpole which pins the flag serving as the reference frame 
flxed to the environment so that the Galilean invariance is violated. However, the 
correct explanation of the flapping flag requires the presence of friction (for recent 
developments see e.g. Ref. jlOOj V 

In the study of coherent quantum systems the special case of a free surface is 
of great interest. Obviously the instability occurs even if the two superfluids are on 
the same side of the interface, i.e. with a free surface over a pool of two or more 
interpenetrating superfluid components, such as a dilute solution of superfluid "^He 
in superfluid '^He or the neutron and proton superfluids in a neutron star (24) . The 
instability exists even in the case of a single superfluid below the free surface under 
the relative flow of the normal and superfluid components, as has been pointed out 
by Korshunov |23| . He also derived two criteria for the instability, depending on the 
coupling to the environment. In this case the frame-flxing parameter which regulates 
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the interaction with the environment is the viscosity 77 of the normal component of 
the Uquid. For ?? 7^ the critical counterflow velocity v — Vs ^ for the onset of the 
surface instability is independent of 77, but does not coincide with the ideal classical 
result which is obtained when viscosity is neglected [rj = 0): 



In laser-cooled rotating atom clouds in the Bose-Einstein-condensed state the insta- 
bility of the free surface is the generic vortex formation process 101, 102 . Another 
case of vortex formation via surface instability is the interface between the normal 
and superfluid states of the same fluid, which are flowing at different velocities. Such 
a situation has been discussed extensively for rotating ■^He-B which is irradiated with 
thermal neutrons ^ I103| . 

Finally, it is worth noting that these ideas might find applications when studying 
the instability of quantum vacuum beyond the event horizon, or the ergoregion of 
the black hole At the superfluid KH instability the ergoregion develops on the 
AB interface, as the energy of the surface waves becomes negative. Such ripplons 
are excitations of the AB interface and provide a connection to the presently popular 
idea in cosmology proposing that the matter in our Universe resides on hypersurfaces 
(membranes or simply branes) in a multidimensional space. Branes can be represented 
by topological defects, such as domain walls or strings, and by interfaces between 
different quantum vacua. In superfluid '^He, the brane is represented by the AB 
interface between the two quantum vacua, 3 He- A and ^He-B. With the AB interface 
instability, it might then be possible to model the instability of the quantum vacuum 
in the brane world. The latter occurs in the ergoregion owing to the interaction 
between the matter on the branes (represented by ripplons) and the matter in higher- 
dimensional space (represented by quasiparticles in the bulk superfluids). 

3.6. Measurement of AB interface instability 

The superfluid KH instability of the AB interface is a new physical phenomenon 
with wide-ranging ramifications, as discussed in Sees. 13.11 - 13.51 It also provides a 
whole new set of tools for further measurements on vortex dynamics. The standard 
KH measurement is that of the critical rotation velocity fie shown in Fig. ^2 Here 
the NMR absorption plotted on the vertical scale monitors the number of rectilinear 
vortices in the B phase while the rotation velocity fl on the horizontal axis is slowly 
increased. Temperature T and barrier current /b are kept constant during the 
measurement. The first discontinuous jump in NMR absorption marks fl^. This 
is the rotation velocity at which the first vortices are formed in the initially vortex- 
free B phase. If is increased further, the instability occurs repeatedly at the same 
value of critical counterflow velocity v(r = i?off) = Vc ~ Vn(-R) — Vsb(^^) = 
(denoted with dashed lines in both panels). Thus the staircase- like pattern of NMR 
absorption is composed of steps of different height and can be calibrated to provide 
the exact number of vortices which is transferred across the AB interface in each single 
instability event. Their number AA^ is denoted next to each instability event. We see 
that AN is a small random number, which after many similar measurements is found 
to be between 3 — 30 with an average of ~ 10 jl[)4| . The measurement in Fig. Illlis 
performed at high temperatures where each of the AA'^ vortex loops rapidly develops 
into a rectilinear B-phase vortex line, after they have been tossed as a tight bundle 
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Figure 11. Measurement of the superfluid KH Instability of the AB interface. 
(Top) A step-like increase in absorption is observed in the Larmor region of the 
NMR line shape when the rotation velocity fl is slowly increased and the critical 
value of counterfiow velocity is reached. The KH instability occurs repeatedly 
in the form of different independent events as a function of fi. A random number 
of AN rectilinear vortex lines are formed in each event in the B-phase section 
of the sample (see Fig. |7J. (Bottom) At the instability the critical counterfiow 
velocity «^c(^off) is abruptly reduced by kAN/ {2itR) to a sub-critical level. The 
constant critical counterfiow velocity Vc gives rise to the dashed lines in each panel. 
iJoff < -R is the effective radius at which the instability occurs (see Refs. llU4ll6bl 
for details.) 



across the AB interface. From such measurements at different temperatures, fl^ has 
been collected to provide the curves at different values of /b seen in Fig. O 

In the measurement of Fig. ^] the sample might be in either the BA or BAB 
configurations. In the setup of Fig.Qthe magnetic field distribution along the vertical 
axis is almost symmetric with respect to the middle of the sample. Also the end plates 
of the sample cylinder and the NMR pick-up coils are at the same distance from the 
AB interfaces in both B-phase sections. As a result the measured values of fie, the 
rotation velocity corresponding to the first instability event, are identical for the top 
and bottom B-phase sections in the BAB configuration. Nevertheless, the events 
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Figure 12. Schematic illustration of how a KH instability event might evolve 
on the AB interface. (Top row) Initially when rotation is increased, vortices 
form at a low critical velocity in the A phase section while none are formed in 
B phase. At the critical B-phase counterflow velocity the AB interface becomes 
unstable towards wave formation, AN vortex loops end up on the B-phase side 
of the interface, and then develop to AN rectilinear vortex lines. The end result 
from the instability is that the A^b = AN vortex lines, which now pass through 
the AB interface after the first instability event at f2c, wind up the B-phase 
superflow velocity and the boundary settles down. In the upper right corner the 
velocities are sketched of the normal (vn) and superfluid {va\, v^b) components 
in the two phases. They are shown here in the laboratory frame as a function 
of r in the situation after the instability event (see also Fig. (Bottom row) 
A schematic illustration of how the vortex injection might happen. When the 
boundary becomes unstable, waves form on the interface. A small number of 
vortex lines becomes trapped in the deepest corrugation which expands to the 
B-phase side (where H < H\b{T) and A phase is unstable). The corrugation 
becomes separated and the A phase shrinks away but the circulation is left behind. 



at the two AB interfaces in the BAB configuration are not correlated and AA^ differs 
evidently randomly from one event to the next. Therefore, as can be seen from Fig. 1111 
subsequent instability events at the two AB interfaces get soon out of step, when 
is increased above fie, and occur in the end randomly, but such that the critical 
counterflow velocity remains fixed: Vc « UcR. The reason for this is that vortices 
play no direct role in the instability condition of Eq. (|3Uf) : Even though the vortices 
in the A-phase vortex layer covering the AB interface are as closely packed next to 
each other as the soft core diameter allows (see Fig. 1131 right), they do not determine 
the instability. Instead it is the stability of the interface itself in the tangential flow 
according to Eq. (|Sn|l . 

A schematic illustration is shown in Fig. 1121 of the process in which AN vortex 
loops might be transferred across the AB interface. When the instability threshold 
is reached while increasing fi, surface waves with wave vector fco = ^F/a begin to 
form on the AB interface. At this stage the AA'' vortices correspond to the A-phase 
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Figure 13. AB interface instability as a function of the restoring force at con- 
stant temperature. (Left) Measured critical rotation velocity Qc of the first KH 
instability event as a function of the current in the barrier magnet. These 
measurements are conducted in the temperature regime of the hydrodynamic 
transition between regular and turbulent vortex dynamics: (a) open symbols 
mark events in which only a small number of AN new vortices are formed in 
the B-phase section; (■) filled symbols denote events in which a turbulent burst 
increases the B-phase vortex number to A^cq. The solid curve represents Eq. lijUl 
without fitting parameters. (Middle) Half of the calculated ripplon wavelength 
A = 2iT/ko at the instability according to Eq. I25i . (Right) The bottom axis 
gives the separation as between vortex quanta on the AB phase boundary on 
the A-phase side, calculated from the number of vortices in the A-phase section 
(when these are distributed evenly in single-quantum units along the outer sample 
circumference r = R). The top axis gives the number of single-quantum units 
expected in one trough (A/2) of the surface wave. This number A/(2as) agrees 
with the measured average for A7V 10 (see Ref. )104l for details). 

circulation quanta which reside in one of the corrugations of the wave where it is about 
to become the deepest and most dominant trough. In the lower row of diagrams in 
Fig. ^1 the possible sequence of events is depicted which might then take place. The 
trough becomes a potential minimum for vortices, which pushes the vortices even 
deeper in this well. Ultimately, as the trough gets distorted, a Magnus force starts 
to act on the vortices. The interface motion is highly damped and it moves faster 
than the vortices. At some point the interface becomes unstable below its equilibrium 
level in the region where H < Hab{T) and it springs back, the trough is smoothed 
out, and a vortex bundle is left behind in B phase. As a result of such a transfer of 
circulation across the AB interface, the counterflow velocity at the phase boundary is 
now subcritical, and the interface settles in its stationary state configuration. 

A further illustration of the KH measurement and its agreement with Eq. 
is shown in Fig. 1131 This plot is generated by measuring Qc in similar manner as in 
Fig. 111! by increasing rotation to fie at different values of /b, keeping temperature 
constant. Alternatively, the measurement can be performed by decreasing /b at 
constant O, so that the critical curve in Fig. El is traversed horizontally from right 
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to left. Since the measurement is at constant temperature, it is also performed 
at constant value of Hab{T, P). The result demonstrates the dependence of the 
instability condition in Eq. H3()|l on the restoring force in Eq. H26|l . owing to the 
dependence of WH on /b. When /b is changed, the AB interface moves within the 
barrier magnet to a location where the field approximates the equilibrium transition 
field, H{r,z) ~ Hab- Because of the large curvature in the solenoidal field and the 
surface tension a of the AB interface, the equality is only approximately followed. At 
large /b the interface is flat, as seen in Fig.d and fie changes only slowly in Fig. El In 
contrast, at low /b the interface is highly curved, its curvature changes rapidly with /b, 
and accordingly Qc is also a strong function of lb- It is interesting to note that the fie 
data with the lowest /b values in Fig. E] correspond already to a configuration where 
the A-phase volume is not a complete layer between two independent B-phase sections, 
but where the AB interface separates out a torus-shaped A-phase volume around the 
cylinder wall, with a B-phase channel connecting the top and bottom sections of the 
sample |1(J5) . Surprisingly, even in this situation the measured data points lie on a 
relatively smooth ilc(-^b) dependence. 

With a given barrier magnet the KH instability is restricted to the range of 
velocities allowed by Eq. . With the barrier magnet of Fig. [3 the critical velocities 
J7c are in the range 0.7 - 1.6rad/s, as seen in Figs. El and El However, vortex 
loop injection into vortex-free B-phase flow can be performed at higher rotation 
velocities using the following procedure: The barrier field is initially reduced below the 
equilibrium value Hab{T, P), so that there is no A phase in the sample volume. The 
rotation velocity is then increased to the desired value above the critical velocity fie 
of KH instability. In this all-B-configuration the vortex free state can be maintained 
up to the velocity at which vortex formation from other sources starts. Next the 
barrier field is ramped up. A phase then forms in a sudden hysteretic transition at 
a magnetic field which somewhat exceeds the equilibrium value Hab- Because of 
this superheated transition, the A phase forms simultaneously within a larger volume, 
which is unstable until the A-phase vortices have been formed. The AB interface is 
also unstable, until a large number of vortices is transferred into the B phase. The AB 
interface finally settles down, when the B-phase superfiow velocity has been reduced, 
usually well below the critical value. This injection technique is useful for studying the 
propagation of a large number of vortices at velocities above fie at any temperature 
where the AB interface exist. 

The KH shear-fiow instability provides a convenient mechanism for precise vortex 
injection into initially vortex-free applied fiow. This is its principal application in 
the dynamical measurements which are described in Sec. 01 Its critical velocity is 
predictable and can be tuned externally. With the KH instability it becomes possible 
to inject a bundle of vortex loops as an externally triggered event: the rotation velocity 
can be stabilized as close as within Afi ~ 0.01 rad/s below the threshold ilc and then 
suddenly increased by 0.02 rad/s to start the KH event. Such reliability in vortex 
injection allows new types of measurements on the dynamic evolution of vortex lines 
in applied fiow. The prime example is the determination of the vortex fiight time, 
a measurement of the velocity at which a vortex propagates into vortex-free applied 
counterflow in a rotating column (Fig.Eti). Even more importantly, the KH instability 
provided the flrst firm identification of the hydrodynamic transition from regular to 
turbulent vortex dynamics which occurs on cooling below ~ 0.6 Tc; in other words, it 
allowed to distinguish a new phenomenon from other sources of uncontrolled vortex 
formation. 
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4. Transition from regular to turbulent dynamics 

4.1. Introduction 

At high temperatures above 0.6 Tc the dynamical behavior of vortices in '^He-B is 
regular, i.e. their number does not increase during a time-dependent process. This 
is seen from the fact that single- vortex steps are observed in Figs. |31 and in the KH 
measurement of Fig. 1111 the number of vortices AA'^ formed in each instability event 
fits within the distribution expected on the basis of Fig. El (right). Some time ago it 
was recognized that at lower temperatures single- vortex processes are not observed in 
large applied flow instead, a large number of vortices is suddenly formed so that 
the final state in rotation appears to include close to the equilibrium number of vortex 
lines. More recently, a consistent explanation of this phenomenon has been presented 

Eg. 

On the basis of extensive measurements on the Kelvin-Hclmholtz instability 
as a function of temperature, it is now understood that mutual friction divides 
the dynamics in '^He-B into a high-temperature regime with regular vortex-number 
conserving motion, and a low-temperature regime where superfluid turbulence 
becomes possible. This hydrodynamic transition is unusually abrupt as a function 
of temperature, as seen in Fig. 1141 This is the reason why it was first mistakenly 
interpreted to signal a strongly temperature-dependent new critical velocity in 
Ref. |3ni- The applied flow velocity, the counterflow velocity Vn — Vg, is an important 
factor in this transition: at low velocity a single vortex has been observed to be in 
stable precessing motion for an entire day at temperatures down to below 0.2 Tc |l()6j . 
A similar result has been verified for ^He-II |1()7| . With increasing flow velocity the 
vortex undergoes an instability towards the formation of Kelvin waves, which may or 
may not lead to a rapid multiplication in the number of vortex lines. The outcome 
from the instability depends on the mutual friction parameter q. 

The transition from regular to turbulent dynamics as a function of mutual-friction 
dissipation is a new phenomenon. It has not been observed in ^He-II; in this case the 
transition is expected only a few tens of fiK below T\, where Ps is vanishingly small 
and the superfluid state very different from that further below T\. In ■^He-B the 
transition is in the middle of the experimentally accessible temperature range where 
the superfluid properties are continuous and well developed. Here the transition can 
be observed in one experiment by scanning temperature from the superconductor-like 
dynamics at high vortex damping to superfluid *He-like turbulence at low damping. 
This shows that superfluid dynamics can be varied and that the traditional "^He- 
like superfluidity is just one limiting case in this spectrum. The opposite extreme 
is superfluid ^He-A where sufficiently low temperatures, to reach turbulence, are 
probably not experimentally realistic. 

In this section we first describe how the transition appears in the KH measure- 
ments. We then proceed with some models for its explanation. 

4.2. Regular vs turbulent dynamics in Kelvin- Helmholtz measurements 

Fig. El shows the KH critical velocity measured as a function of temperature. Here 
the data points have been classified according to whether the final state after the 
first critical event only includes the vortices generated in the KH event itself (open 
symbols), or if it incorporates almost the equilibrium number of vortices (filled 
symbols). A sharp division line, with little overlap of open and filled data points. 
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Figure 14. Critical rotation velocity Qc of the first KH instability event at 
high barrier field, when the A-phase section is stable down to T — > 0, or 
-ffb(''i •^)|max > Has{T = 0). The data points have been classified according 
to the nature of the dynamic evolution after injection: (o), regular vortex number 
conserving; {•), turbulent burst. A transition in the dynamics is seen to take place 
at about 0.59 Tc, with little overlap of open and filled symbols in the transition 
regime. The solid curve denotes the calculated dependence from Eq. l3Ui without 
fitting parameters. In the absence of more appropriate parameter values, mainly 
isotropic zero-field values have been used which do not seem to provide perfect 
agreement at high fields and low temperatures. 



is seen to divide the plot into two regimes as a function of temperature: at high 
temperatures only the KH vortices are generated while at low temperatures the 
equilibrium vortex state is obtained after the instability event. The KH critical velocity 
as a function of temperature is a continuous smooth curve across this division line and 
continues to follow the calculated dependence. The same features are illustrated in 
Fig. 1131 which is measured at a temperature close to the division line. A surprising 
and characteristic property of these plots is that there are essentially no data points 
with an arbitrary intermediate number of vortices, not even in the transition regime. 

The conclusion from these measurements is that at temperatures below the 
division line a short turbulent burst follows the KH event, after the closely packed 
bundle of vortex loops has arrived across the AB interface into the vortex- free B-phase 
flow fFig. I12|l. The turbulent burst generates the vortices needed to reach almost the 
equilibrium number of vortices for the B-phase flow. The propagation of these vortices 
along the long rotating sample as a vortex front followed by a twisted vortex bundle 
is discussed in Sec. |5| A compelling argument for the interpretation in terms of a 
turbulent burst is the continuity of the KH critical velocity across the division line - 
it is unrealistic to assume that the nature of the instability would change so suddenly. 
Since these early observations it has been understood that the KH instability 
is just one example of vortex seed loop injection in externally applied flow. Other 
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Figure 15. The number of rectilinear vortex lines N after KH injection, 
normalized to the equilibrium number of lines A^eq, plotted as a function of 
temperature. At around 0.58 Tc, a sharp transition in the number of lines is 
observed. At high temperatures the injection results in only a few lines, but at 
low temperatures the equilibrium number of lines is observed 11041 . 



examples are examined in Sec. 14.61 However, since the KH instability was originally 
the most convincing case of vortex injection and of the turbulent burst, with very 
particular properties, we describe this example here in more detail. 

Figure [Tsl shows a close-up of the transition region as a function of temperature. 
The number of vortex lines N is plotted, normalized by the equilibrium number A'eq, 
after a single event of KH injection. At high temperatures the injection leads always 
to a small number of rectilinear lines, as expected for the KH instability fFig. 
Below 0.6 Tc the number of lines, after the system has settled down, is very different: 
the injection results in the almost complete removal of counterflow. The two regimes 
are separated by an abrupt transition at Ton ~ 0.58 Tc which has a narrow width 
of ATon « 0.04 Tc. We attribute the final state at temperatures below Ton as the 
fingerprint from a brief burst of turbulence which multiplies the vortex number close 
to iVeq. Here with KH injection the turbulent burst occurs at the injection site at a 
distance of 30 mm from the closest end of a detection coil. Thus there is no direct 
NMR signal which would identify the burst itself, only the propagating vortex state 
after it has traveled from the AB interface to the detection coil. This means that the 
turbulent burst is short in duration and localized to some section of sample length 
adjacent to the AB interface. 

In the transition regime the turbulent burst may not be triggered by the first 
vortex injection event. Instead, a turbulent burst can be preceded by one or more KH 
injection events which do not lead to vortex multiplication. However, once turbulence 
sets in, it generates the equilibrium number of lines in almost all the cases. Such 
events where turbulent vortex multiplication terminates before the number of lines 
reaches the equilibrium value are rare. 
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Figure 16. Temperature-velocity pliase diagram of turbulence. Each marker 
in the plot represents the result from a measurement where rotation is increased 
from zero to the critical velocity Qc of the AB interface instability. The cases 
where turbulence is observed after the instability are marked with filled symbols 
(I) and regular cases with open symbols (□) ^25j. The horizontal top axis 
gives the temperature dependent and velocity independent dynamic parameter 
q = a /(I — a'). We conclude that the phase boundary is primarily determined by 
its mutual friction dependence. 



The transition from Fig. ^] is shown as a phase diagram in Fig. ^1 plotted 
in terms of the B-phase counterflow velocity v — flR and temperature T. In this 
diagram each data point represents a KH injection measurement, accumulated with 
different settings of the externally controlled parameters, so that as wide a variation 
as possible is obtained for the critical rotation velocity flc and temperature T. Each 
marker in Fig. 1161 thus indicates a KH injection event into vortex-free counterflow 
with some parameter values which are not of interest in this context. What we are 
interested in here is the division in filled and open symbols: events followed by a 
turbulent burst are again marked with filled symbols (■) while events which only lead 
to a few rectilinear lines are marked with open symbols (□). The transition occurs at 
about 0.59 T^, at higher temperatures no injection events lead to vortex multiplication 
while at lower temperatures all injections lead to the equilibrium vortex state. The 
striking conclusion from this plot is that the phase boundary between turbulence at 
low temperatures and regular dynamics at high temperatures is vertical and thus 
independent of the counterflow velocity. 

Plots similar to Fig. 1161 were also measured at 34 and 10 bar pressures |108) . 
In Fig. El the results are summarized. We assume that the transition is velocity 
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Figure 17. Transition between regular and turbulent vortex dynamics at different 
pressures. The transition is assumed velocity independent and the events are 
categorized according to their temperature I1U8I . The transition temperature 
increases with pressure, but the half width of its distribution is approximately 
<j ~ 0.03 Tc at all pressures. 



independent in the measured range of velocities, and compile the data to show the 
probability of the transition between regular and turbulent cases as a function of 
temperature. The fitted Gaussian distributions give a narrow half-width of cr « 0.03 Tc 
at all pressures, centered around a transition temperature in the range 0.52 - 0.59 Tc, 
depending on pressure. Mutual friction data is available at 10 and 29 bar pressures 
|72|. At lobar the transition occurs at q = 0.6 and at 29 bar q — 1.3. Thus the 
transition appears to move to higher q value with increasing pressure. Measurements 
were also carried out at zero pressure where the transition was found to be below 
0.45 Tc at low rotation velocities of 0.5 - 0.7rad/s. 

To conclude, KH injection measurements indicate that the phase boundary 
between turbulent and regular vortex dynamics is foremost a function of temperature 
and independent of the applied flow velocity at velocities above |v| = |vn — Vs| > 
2.5mm/s. Thus the transition occurs as a function of mutual friction, such that the 
dynamic parameter q is of order unity. 



4-3. Classical and superfluid turbulence 

Since classical and quantum turbulence share many common features, we begin with 
the basic concepts from classical turbulence by inspecting the properties of the different 
terms in the Navier-Stokes equation |M] 

— + V/i = vxu; + j/V^v . (37) 
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Here v is the kinematic viscosity (viscosity 77/density p), u; = V x v is the vorticity in 
the inertial (laboratory) reference frame. Turbulence is governed by the interplay of 
the two terms on the r.h.s. of this equation, the inertial (first) term and the viscous 
(second) term. 

The transition to turbulence is determined by the Reynolds number Re = LU/u, 
formed from the characteristic values for the three quantities describing the flow: its 
velocity J7, the geometric size of the system L, and the kinematic viscosity v. For small 
Reynolds numbers, the viscous term, — i/fc^Vk for a perturbation with wave vector k, 
stabilizes laminar flow. In contrast, at large Re 3> 1 the inertial term in Eq. (|37|l 
dominates, and laminar flow becomes increasingly unstable against the formation of a 
disorganized flow of eddies. In the most carefully prepared experiments laminar flow 
has been maintained in a circular pipe up to Re ~ 10^ However, the higher the 

Reynolds number the smaller the disturbance needed to trigger turbulence IllOj . The 
evolution of turbulence is governed by the Kolmogorov energy cascade: the kinetic 
energy of the flow is transferred to smaller and smaller length scales, via the decay 
and break up into ever smaller vortex loops along the so-called Richardson cascade 
until a length scale is reached where the energy can be dissipated by viscosity. 

In superfluids turbulence acquires new features. First, the superfluid consists of 
two inter-penetrating components: the frictionless superfluid and the viscous normal 
fractions. The total density p ^ pn + Ps is the sum of the densities of these two 
components. The normal component behaves like a regular viscous fluid while the 
superfluid component is an ideal superfluid. At the superfluid transition the density of 
the superfluid component vanishes, but increases with decreasing temperature, until 
in the T — > limit the normal component vanishes. Secondly, the vorticity of the 
superfluid component is quantized as discussed in Sec. 12.11 If both the normal and 
superfluid components are able to move, the turbulent state bears more resemblance 
to the turbulence of classical viscous liquids. This is often the case in superfluid ^He-II 
In superfluid "^He, however, the normal component is so viscous that it is essentially 
immobile. The flow is then carried by the superfluid component which contains a large 
number of quantized vortices. Here a new class of turbulent flow becomes possible: 
one- component superfluid turbulence. It is this state of turbulent flow that we consider 
in what follows. A window for this regime of superfluid turbulence appears because 
of the existence of one more large parameter in superfluid ^He, namely the ratio of 
its normal kinematic viscosity to the superfluid circulation quantum, v / k. In '^He at 
T ^ 0.5 Tc this ratio is v/n ^ 10^ in contrast to ''He-II where it is of the order of 
unity. A more detailed discussion is given later in this section. 

Because of its large viscosity, the normal component of '^He-B moves together 
with the container. As a result, Eq. (|ll|l takes the following form in the frame where 
the normal component is locally at rest (v,i — 0): 



is the viscous part of the mutual friction force. 

The flrst inertial term on the r.h.s. drives the flow instability towards turbulence 
in the same way as the inertial term in the Navier-Stokes equation does for potential 
flow in classical hydrodynamics (HH. The second dissipative term on the r.h.s. of 



dt 



+ 



V/i = (1 - a')vs X + fvi 



(38) 



where 




Dynamics of vortices and interfaces in superBuid ^He 



44 



Eq. (|38|l is the counterpart of the dissipative term in the Navier-Stokes equation H37|l. 
It tends to stabilize the flow, since it leads to energy dissipation 

d Vg 

where Vg • fvisc < according to Eq. H13|l . The fundamental difference between the 
dissipative terms in classical and superfluid dynamics is that in superfluids this term 
has the same scaling dependence on velocity and vorticity /vise ~ ctoJsVg as the inertial 
term. This is a consequence of the two-fluid dynamics, where the vortices provide the 
mechanism of momentum and energy transfer between the two components of the 
fluid. Thus the effective Reynolds number - defined as the ratio of the inertial and 
dissipative terms in the relevant hydrodynamic equation - has to be changed. In 
superfluids, it is the ratio of the two terms on the r.h.s. of Eq. (|38|l : 

= (1 = 1/g ■ (40) 

According to this definition, and in analogy with classical hydrodynamics, turbulence 
in superfluids is expected when !/(? ^ 1 and laminar (regular) flow when 1/q < 1. 
As distinct from viscous liquids, this condition is independent of extrinsic quantities 
such as the counterflow velocity or the characteristic dimension R of the sample. The 
temperature dependence of 3? for '^He-B and He-II is explicitly shown in Fig. [S] 

The above dimensional arguments on the transition between regular and turbulent 
vortex dynamics are in agreement with the results in Fig. 1161 where the boundary 
between the laminar and turbulent regimes was also found to be at g ~ 1. Thus the 
velocity independent parameter 3? = 1/q, which controls the transition to turbulent 
flow in superfluids, plays the same role as the velocity dependent Reynolds number 
Re = UR/v in classical hydrodynamics. 

These considerations are valid provided that the applied counterflow velocity is 
high enough to sustain vortices, i.e. it exceeds the Feynman critical velocity fSec. l2.3|l . 
which is easily satisfied in any experiment. From the Feynman criterion, one can define 
another dimensionless parameter, the "superfluid Reynolds number" 

RCs = U^R/n, (41) 

where Us is the mean superfluid velocity with respect to the normal component, i.e. 
the counterflow velocity. If the condition Rcs > 1 is fulfilled up to a logarithmic 
prefactor ln(i?/^), it becomes energetically favorable to add a vortex line in the bulk 
flow. 



Onset of turbulent hurst 

In KH injection a tight bundle of many small vortex loops (Fig. I21|) is transferred 
across the AB interface into vortex-free B-phase flow. How is a turbulent burst started 
from such a vortex bundle? In this section we discuss the initial phase of vortex 
multiplication in applied flow using a simple phenomenological model from Ref. jll2| 
which is constructed in the spirit of the Vinen equation for superfluid turbulence |118j . 

To characterize the initial conditions we need two numbers: the intrinsic velocity- 
independent 3? = 1/q in Eq. (|in|) and the "superfluid Reynolds number" Res in 
Eq. |0lj- We assume Rcs ^ 1, which corresponds to the typical experimental situation 
in ^He-B that the presence of many more vortices would be energetically possible. If 
a large energy barrier prevents vortex nucleation, then vortices are not necessarily 
created even at high Res ^ 1 and the superfluid remains in a metastable state of 
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counterflow. At velocities well below the intrinsic critical velocity of vortex formation 
(Sec. 12 .^11 . Vc ~ k/£, superfluid turbulence can be initiated if quantized vortices 
are injected by some extrinsic means into vortex-free flow. In the rotating sample, 
provided that 1/q > 1, turbulence then develops in a sudden burst where the initial 
seed vortices start to multiply and form a vortex tangle in the bulk volume. 

We start from an initial configuration containing many (essentially more than one) 
randomly oriented vortex loops, which might have been either injected or created by 
some other precursory process from a set of a few seed vortices. An example of the 
precursor can be the Kelvin-wave instability discussed in Sec. 12.61 We show that 
the initial array of entangled vortices is unstable towards a burst-like multiplication 
of vorticity provided the mutual friction is low enough. We assume that vortex 
multiplication occurs within a certain region in the fluid where the number of vortices 
is large. The region where multiplication takes place can exist near the location of 
the seed loops and/or near the walls of the rotating container, where the rotating 
counterflow velocity reaches its maximum value and the number of vortices in effect 
is increased by the presence of image vortices. The vortex tangle created in such a 
region of multiplication penetrates next into the rest of the fluid volume. 

Since superfluid vorticity is quantized, the formation of new vortices during the 
onset becomes the key issue. We consider this process, taking mutual friction into 
account, and derive an equation for the evolution of the density of entangled vortex 
loops during the initial stages of the transition to turbulence. The multiplication of 
seed vortices can be studied with Eq. . Taking the curl of both sides of Eq. 
we obtain the equation for the superfluid vorticity 

^ = (1 - a') [V X (v, X a;,)] + V x f^isc • (42) 

Here we assume that the vorticity produced by turbulence is much larger than the 
angular velocity, ujs S> f^, i.e. the vortex density is much higher than that in 
equilibrium. This allows us to neglect V x Vn = 2f2 in comparison to Wg. 

Let us now average Eq. over randomly oriented vortex loops with dimensions 
spread over an interval around a characteristic size i, using the ideas of Ref. |113| . 
In a state of entangled vortex loops, their three-dimensional density is Uv 
while the vortex-loop length per unit volume (two-dimensional vortex density) is 
L = £n^ = = nv^^ After averaging Eq. (|42|l only its scalar value is meaningful, 
as any of one its components to a specific direction vanishes. Let us express the two 
terms on the r.h.s. of Eq. (|42|l in terms of the vortex density L, keeping in mind 
that the velocity produced by the vortex tangle is of order k/£, while its vorticity is 
Us ^ a/i"^ ~ HL, where k = (k/Att) \n{i/^). 

Regarding the onset of turbulence, the reactive coefficient a' in Eq. (|42|l simply 
renormalizes the inertial term of conventional hydrodynamics [the first term on the 
r.h.s. of Eq. 1(23] that drives the instability towards turbulent vortex formation. 
Therefore, the vortex density increases owing to the first term in Eq. I)42|l according 
to 

L+ = A{1 ~ a')vsL^/'^ = A{1 - a'){U - vq)L^''^ (43) 

where ^ ^ 1 is a constant. The superfluid velocity in Eq. (|43|l is assumed to be 
Vs — U ^ V{) where U is the counterflow velocity, and ~ k,/£ is the self-induced 
velocity for a vortex loop of length t and core radius ^. The kinetic energy of the 
superfluid grows due to the increase in the loop density. The energy is taken from the 
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external source at the length scale R with the rate dE/dt ^ Ei^L^, where Ei^ is the 
energy of the vortex per unit length. 

Vortex multiplication saturates when the loop density reaches a value such that 
U — vq kL^/'^. If the density happens to become larger, it will decrease towards 
Lsat ^ {U /hY' while the kinetic energy is returned back to the external source. In 
other words, saturation is reached when the "turbulent superfluid Reynolds number" 
j^g(tuib) [/^^^ becomes of order unity. The condition Rcg S> 1 ensures the separation 
of scales, £ ^ R, that is required for the formation of a vortex tangle. The vorticity 
at saturation, tJsat ^ '^/^sat ~ Rcg, is much larger than the equilibrium vorticity 2fl. 

The second (dissipative) term in Eq. (|42|l acts to stabilize vortex-free flow, thus 
reducing vortex density in a way similar to that in normal fluids. Estimating the 
dissipative term as /vise ~ ctujsVs we find for the rate of decrease in vortex density 

L.^-Ba{U -vo)L^''^ (44) 

where i3 '--^ 1 is a constant. In contrast to viscous fluids, the dissipative term has here 
the same scaling dependence as the driving term in Eq (|43|l . 

The multiplication of vortex loops, as described by Eq. H43|l . can be understood 
in terms of vortex collisions and interconnections. Such processes were indeed seen in 
numerical simulations on quantized vortices [251 11141 180| . Reconnections of vortices 
accompanied by the formation of a vortex tangle in normal fluids were considered 
recently in Refs. jll5l Each reconnection of quantized vortices takes place 

over a microscopic time of the order of the quasiparticle collision time, which is 
much shorter than the characteristic times involved in hydrodynamic processes. It 
is accompanied by some small amount of dissipation within a volume of the order 
of We consider these reconnection processes as instantaneous and neglect the 
corresponding dissipation. The rate of increase in vortex loop density should be 
quadratic in riv, and thus ri+ = Av-cU^l"^ . Here Ur is the relative velocity of the vortex 
loops, £^ is the loop "cross section" , and the constant A ^ 1 describes the "efficiency" 
of the vortex multiplication due to pair collisions. Using the definition of L the vortex 
multiplication rate becomes L+ ~ v-^L^/'^. The vortex velocity is determined through 
the mutual friction parameters a and a' such that vl = (1 — ol')ws — atl^s x Vg. After 
averaging over randomly oriented vortex loops the last term vanishes, resulting in the 
average relative velocity of loops Wr proportional to the longitudinal component of vl , 
Vj: (1 — a')vs. The rate L+ thus agrees with Eq. H43|) . 

The effect of the (second) viscous term on the r.h.s. of Eq. (|42|l is to decrease 
the loop density in the region of multiplication, by inflating the loops in the applied 
counterflow, and by extracting them from the region of multiplication into the bulk 
where the counterflow is smaller. The viscous component of the mutual-friction force 
leads to variations in the vortex loop length £ ^ 27rwL ^ otVs ^ a(U — vq). Finally, for 
the rate of variation in the vortex-loop density from the viscous component one obtains 
Eq. (|44(l . The length increases while the density decreases as long as saturation is not 
reached, U — vq > 0. If saturation is exceeded, vq > U, the density increases since 
the loops shrink owing to the friction. Accordingly, the vortex loops are extracted 
from the region of multiplication into the bulk before saturation U > vq, and they are 
extracted out of the bulk fluid \iU < vq. 

As we see, both inertial and viscous mutual friction terms, Eqs. (|43ll and 144() . 
have the same dependence on vortex density, i.e. on the vortex length scale. The 
total variation in loop density in the region of multiplication is the sum of the two 
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processes, L = i+ + L_. Putting vq = k/i = kL^^^ we obtain 

{l-aT'L = iq,-q)(uL^^^~kL^) . (45) 

Here the critical value Qc — A/B is generally of order unity. 

Equation (|45|l resembles the Vinen equation 113 for turbulence in superfluid ^He- 
ll. However, the difference is that the coefficient qc — q can now have either positive 
or negative sign, depending on the mutual-friction parameters. As a result, two limits 
can be distinguished. 

In the low-viscosity regime, which is typical for "^He-H, q^ — q> Q. In this regime, 
the rate of multiplication is faster, and the number of created vortex loops is large: 
each new vortex loop serves as a source for producing more vortices. As a result, 
avalanche-like multiplication takes place, which leads to the formation of a turbulent 
vortex tangle. As the number of vortex loops grows, the self-induced velocity increases 
and finally the saturated density Lgat is reached. 

In •^He-B the opposite regime is possible, with qc — q <Q. In this viscosity-domi- 
nated regime the rate of extraction of vortex loops exceeds the rate of multiplication; 
there is no time for vortices to multiply since all seed vortices and newly created 
vortices are immediately wiped away into the bulk fluid. Provided initially created 
vortex loops have the density L smaller than the saturation value Lgat ~ (U/k)'^, 
the number of vortices in the final state is essentially equal to the number of initial 
vortices, and the turbulent state is not formed. The corresponding stable solution to 
Eq. H45|l is L — > 0. Equation (|45|l in this limit describes also one more regime of the 
vortex tangle evolution which is realized when the initially created vortex loops have 
the density L, exceeding Lgat ~ (U/k)'^: The vortex loops will shrink down to very 
small sizes due to the viscous mutual- friction force and collapse, £ while L ^ oo, 
since the counterflow U is no longer sufficient to support loops of such small size. In 
this regime, the solution of Eq. H45(l approaches another point of attraction, L — )■ oo. 

One can see that the condition of instability for an entangled vortex array 9 ^ 1 is 
essentially the same as the condition for the propagation of underdamped Kelvin waves 
along an isolated vortex line, established in Sec. 12.61 This indicates that the threshold 
g < 1 in the mutual-friction controlled dynamics is not just an accidental combination 
of parameter values, but may be of more general and fundamental importance for 
superfluid dynamics; however, its role remains yet to be fully investigated. 

To justify our assumption that the normal component does not participate in 
superfluid turbulence we compare the viscous force rj^k^Vn in the Navier-Stokes 
equation H37(l and the mutual friction force of Eq. (|12|l exerted on the normal 
component _Fmt Q!psiJ-'s(i^s — v^)- Here rj^ — p^iy is the normal dynamic viscosity, 
and k is the wave vector of the velocity field. Variations of Vn are smaller than those 
of Vs when vk^ ^ [p^/ Pn)aoJs- We estimate k^ ~ and ^ k/I"^ in terms of the 
vortex-line density. The condition becomes 

i^/k > {ps/pn)a ■ 

On the left-hand side of this inequality we find a new parameter, the ratio of the 
kinematic viscosity and the circulation quantum v/n, which is a characteristic of 
the particular superfluid. For "^He at T = 0.5 Tc we have v = rj^/ p^ ~ 0.4cm^/s 
and K=:0.066 mm^/s, so that v/k ^ 10'^. The inequality is then well satisfied since 
Ps/pn ^ 1 and a ~ 0.5 for T = 0.5 Tc. Therefore, the normal component remains 
at rest in the container frame. Note that in ^He the situation can be qualitatively 
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different, since the normal component can be involved in the fluid flow owing to its 
much smaller viscosity. 

The overall evolution of the vortex density can be seen as an interplay of two 
processes. The first is the turbulent instability in the region of multiplication, governed 
by Eq. (|43|l . The second process is the expansion of vortex loops from the region of 
multiplication into the bulk due to the dissipative component of the mutual friction 
force. Eq. I|44() taken with the opposite sign gives thus the rate of vortex-loop-density 
flow into the bulk. In this form, it exactly coincides with the Vinen equation as 
derived by Schwarz |114| (see also IllTp . and applies to bulk superfluid turbulence 
which is continuously sustained by an external source that acts in the region of 
multiplication. A typical example is the grid turbulence considered in Refs. |113llll4j : 
the region of multiplication is here assumed to be the vicinity of the grid. Therefore, 
the rate at which vortices are supplied into the bulk from the region where they are 
generated is 

ibuik = +Ba{U - vo)L^^^ = Ba{UL^''^ - kL'^) . (46) 

If the supply continues long enough, the vortex density in the bulk also saturates at 
Lsat = {U/k)^. The solution |118] describing the relaxation towards saturation for 
constant U has the characteristic rate r^^ — aU"^ /2k,. This relaxation time increases 
with decreasing temperature as the parameter a decreases. 



^.5. Energy cascades in developed superfluid turbulence 

In the previous section 14.41 we discussed the onset of superfluid turbulence and 
the resulting state which is characterized by a single scale £ = k/U at which 
saturation occurs. We call such a single-parameter state the Vinen regime of superfluid 
turbulence. This state is very different from the turbulence in viscous liquids, where 
the energy spectrum obeys the celebrated Kolmogorov-Obukhov 5/3-law 

Ek - Ce'/^fc-5/3 . (47) 

Here E^ is the one-dimensional density of turbulent kinetic energy in fc-space, defined 
such that the total energy density E (in physical space) is given by 

£;= l(|v|2) = j dkEk, (48) 

and £fc is the energy flux in /c-space, which is constant in the inertial range of k, where 
viscous dissipation can be neglected, = e. 

Here, we study what is the outcome of the turbulent instability in a superfluid 
in different regimes of Reynolds parameters ^ — \/q and Re^ . We show that in the 
intermediate range of Reynolds numbers 1 /q (Fig. I18|) the turbulent instability leads 
to a state of developed turbulence which is closer to its classical analogue: It exhibits 
a Richardson-like cascade and is different from the Vinen state of turbulence. 

As a starting point we utilize the coarse-grained hydrodynamic equation for the 
dynamics of the superfluid with distributed vortices, Eq. (|16|l . or, after taking the curl, 

^ = (1 - a')V X (v, X u,) + aV X x (a;, x v,)] . (49) 

As we have seen, turbulence develops only if friction is relatively small compared to 
the inertial term, 1/q ^ 1. Here, we discuss the regime of well developed turbulence 
which occurs at5R = l/g^l when the inertial term is strongly dominating. In this 
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limit, g « a, while both a' <C 1 and a <C 1. We show that well developed turbulence 
can have an analog of the Richardson-Kolmogorov cascade, which becomes modified 
by the nonlinear mutual-friction dissipation. 

The main difference of this cascade in superfluids from that in viscous liquids is 
that the dissipation of energy from the mutual friction in Eq. (|39|l occurs at all scales, 
and thus the energy flux Sk must be essentially /c-dependent. From Eq. 139(1 it follows 
that the energy losses from dissipation are 

-^=-VEk , r^g^o , (50) 

where wq = is the average vorticity. In steady-state turbulence, these energy 

losses must be compensated by the energy exchange between different k in the cascade, 
dek/dk = —TEk- Using the relation between the energy Ek and the energy flux in 
momentum space in Eq. H47|l , which follows from general dimensional reasoning in 
the spirit of Kolmogorov, one obtains the following balance equation for the flux Sk' 

^ = -r4/Vv^ (51) 

Such an equation for the energy budget has been used in Refs. |12U[ I12ll I122| : the 
more complicated version with a second derivative |123[ I124| was used for superfluid 
turbulence by Vinen ^25 . In the absence of dissipation, i.e. when F = 0, Eq. (|51|l 
immediately produces the solution with constant energy flux, Sk ^ s. Then Eq. (|47|l 
turns into the Kolmogorov-Obukhov 5/3-law for E^ ~ e^/'^fc"^/'^. 

Eq. 1)5 l|l must be supplemented by a boundary condition: a fixed energy influx 
into the turbulent system from large length scales of order of container size, e^^i/jj = 
U^/R, where U is the counterflow velocity at this scale. A general solution of Eq. ((51|1 
gives the following energy spectrum: 



1+ Jn9/. -7 



(52) 



k{kRy/3 [ {kR)V3 
where the dimensionless parameter 

7 = TR/U = qujoR/U , (53) 
and the mean vorticity is expressed through the energy spectrum, 

ujl = {\uj\f^ f dkk'Ek. (54) 



The ultraviolet cut-off /c* in Eq. H54|l is determined by the microscopic scale at which 
the circulation in the fc*-eddy reaches the circulation quantum k: Wfc^/fc* — k, and 
thus the coarse-grained dynamics is no longer applicable. Since k'^vf. — Ek, the cut-off 
is determined by the spectrum and, as a result, we obtain a closed system of equations 
which can be analyzed for different regions of the Reynolds parameters, = l/g and 
Re, = RU/n. 

Let us consider the turbulent state that corresponds to g <C 1 and g^Re, ^ 1. In 
this case the parameter 7 is close to unity, 1 — 7 ~ (g^Reg)"""^/^ <C 1, and the solution 
H52|l has the form 
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— ■■ ► 

q = qc 



Figure 18. Possible phase diagram of dynamical vortex states in the (g, Res) 
plane. At large flow velocity Reg 3> 1 well above the Feynman critical velocity, 
the boundary between turbulent and 'regular' vortex flow approaches a vertical 
line q = Qc ^ 1. The dashed curve marks the crossover between two regimes of 
supcrfluid turbulence occurring at small q: (i) Developed turbulence of classical 
type, which is characterized by two Richardson-type cascades owing to mutual 
friction dissipation. The Kolmogorov-Obukhov law oc k~^^^ coexists with the 
Ek oc law. (ii) Quantum turbulence of the Vinen type at even smaller q, 
which is characterized by a single length scale £ = k/U. 



where fccr marks the crossover between the Kolmogorov law ai k > ka to the steeper 
law Ek oc k~^ at k < k„. Three important scales - the scale k ^ \/R at which 
pumping occurs, the crossover scale fccr, and the microscopic (quantum) cut-off scale 
/c* - are well separated in this regime: 

fc* - fccr<Z"^/' » fccr ~ 4(Re,(72)-l/2 » }_ . (55) 
ii ri 

At g^RCs ^ 1 one has — 1/R, i.e. the region of the fc~^ spectrum shrinks. 
Here two scenarios are possible. In the first mutual friction is unessential and thus is 
unable to compensate the Kolmogorov cascade. When the intervortex distance scale 
is reached the Kolmogorov energy cascade is then transformed to the Kelvin-wave 
cascade |5] of isolated vortices. In the second scenario the turbulent state is completely 
reconstructed and the Vinen state discussed in Sec. 14.41 and Refs. |113II114| emerges. 
This state contains a single scale £ = k/U and thus no cascade. A possible phase 
diagram of the turbulent states is shown in Fig. ^1 The connection of this phase 
diagram with the flow states observed in various experiments on superfluid ^He-II and 
^He-B is discussed in Ref. |T^ . 

These phenomena found in ^He-B added a new twist to the general theory of 
turbulence in superfluids which was developed earlier by Vinen |127l |2] and which 
was based on numerous experiments in superfluid ^He-II where the first signs of 
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Figure 19. Principle of vortex injection measurements. Different techniques have 
been explored to inject vortex seed loops in rotating vortex- free flow. One can 
then study the evolution of the seeds as a function of temperature. Two regimes 
are found; (1) at high temperatures the vortex number stays constant while (2) at 
low temperatures a turbulent burst increases the vortex number close to that in 
the equilibrium vortex state. The exact onset temperature Ton of turbulence has 
been found to be sensitive to the details of the injection process. This indicates 
that the energy barriers of the different processes, which destabilize the seed loops 
and lead to their turbulent multiplication, depend in addition to temperature also 
on the number, configuration, and density of the loops. 



turbulence were observed already in the 1950's jllHj . The new theory based on ■^He- 
B experiments incorporates two Reynolds parameters (the velocity-dependent UR/k, 
and the velocity-independent q). It suggests different types of developed superfluid 
turbulence in different regions of Reynolds parameters, and allows to derive deviations 
from the classical Kolmogorov-Obukhov scaling law E{k) ~ k^^/'^ . 

The extension of these ideas to the more general case when both the normal and 
the superfluid components may become turbulent has been performed in Ref. |128| 
on the basis of two- fluid hydrodynamics. The results of this analysis are applicable 
to superfluid ^He-II, where the viscosity of the normal component is many orders 
of magnitude smaller than in '^He-B, or in ^He-'*He mixtures where, owing to the 
presence of '^He quasiparticles, mutual friction at low temperature is also expected to 
be significantly higher than in pure superfluid *He-II. 
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4-6. Injection of seed vortex loops in applied counterflow 

Why do seed vortices injected in applied counterflow undergo a transition to 
turbulence? A superficial answer can be given by inspecting the diagram in Fig. [T^ 
Here vortex-free flow is generated by subjecting the superfluid sample to uniform 
rotation. This is a metastable state of high energy where the superfluid fraction 
does not participate in the rotation. As long as the resulting maximum counterflow 
velocity remains below the limit of spontaneous vortex formation of the particular 
sample setup, the vortex-free state persists and no vortices are formed. By injecting 
vortex seed loops in the applied flow the energy barrier for the nucleation of vortices 
can be externally bypassed. We can then watch what happens to the injected seeds. 
At high temperatures damping is large, the seed vortices quickly evolve to rectilinear 
lines, the number of vortices remains conserved, and the flnal state remains highly 
metastable. 

However, below a sudden onset temperature Ton the situation changes radically: 
dissipation in vortex motion has dropped to a level where Kelvin wave excitations 
are not over-damped and helical waves with large amplitude can be formed in such 
sections of the seed vortex where the flow has a component oriented parallel to the 
vortex core. The formation of new vortices in dynamic processes becomes then 
possible and suddenly the superfluid can reduce its overall energy state and reach 
equilibrium. The bottleneck in this chain of instabilities is the single- vortex instability 
in applied counterflow. However, eventually at sufficiently low temperature even this 
will inevitably happen. If the density of the injected seed loops is high enough from the 
beginning, as appears to be the case in KH injection, then inter-vortex interactions 
allow turbulence to switch on immediately and at a higher temperature than the 
single-vortex instability. The evolving turbulence becomes a collective phenomenon, 
which can be described either in terms of the Vinen-type theory in Sec. 14.41 or the 
developing superfluid turbulence in Sec. 14.51 Unfortunately, in rotating flow of ^He-B 
above 0.4 Tc turbulence is limited to a very short burst so that a measurement of its 
evolution with time becomes a challenging effort. The propagation of the vorticity in 
the rotating cylinder after the turbulent burst is described in Sec. [3 

Measurements employing injection of seed loops in applied flow are a new feature 
in turbulence studies. They have become possible in ■^He-B where vortex-free flow 
can be achieved at sufficiently high flow velocity. With the aid of seed loop injection 
one can study different instabilities in the serial chain of dynamic processes which 
lead to vortex multiplication and ultimately to turbulence. One may expect that 
the closer the injected seed loop conflguration is to a situation where the loops can 
interact and create a turbulent burst the higher will be the transition temperature 
to turbulence. Injections with a small number of loops which are far apart require 
a lower temperature to develop to turbulence, i.e. the transition becomes more and 
more irregular being dependent on such details as the number of injected loops, their 
shape, size and density, their mutual interaction and interaction with the container 
walls, where also surface roughness and pinning may matter. These dependences can 
be investigated using different injection techniques |129j . 

Magnetically driven B—fA transition: Large numbers of seed vortices can be 
created by sweeping up the barrier fleld in the setup of Fig. [7| The fleld sweep 
is conducted in rotation at constant temperature, when the entire sample is in B 
phase and vortex free. When the fleld reaches some slightly hysteretic value above 
Hab{T,P) over a short section in the middle of the long sample, then the A phase. 
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its vortices, and two AB interfaces all start to form essentially simultaneously jl()5j . 
A detailed explanation how this happens (in a situation which is unstable from the 
beginning) has not been worked out, but if the rotation velocity f2 significantly exceeds 
the critical KH velocity flc(T,P) a large number of vortices is suddenly injected in 
the B-phase sections. 

Kelvin- Helmholtz instability: The Kelvin-Helmholtz instability of the AB inter- 
face, discussed in Sees. 13.61 and 14.21 is so far the most reproducible of the available 
vortex-injection mechanisms. The injection is carried out at constant external con- 
ditions while the rotation velocity O is suddenly incrementally increased by Af2 « 
0.05rad/s across the critical KH value nc{T,P). At temperatures above Ton such 
injection generates a limited number of vortices in the B phase and consequently the 
staircase pattern in Fig. can be displayed, by triggering multiple instability events 
one after the other. Below Ton the staircase pattern is not obtained, since already 
the very first injection sends the sample in the equilibrium vortex state and removes 
essentially all applied counterflow. The measured signature from the injection in the 
turbulent temperature regime is shown in Fig. 1201 The radical difference from Fig. 1111 
is evident: The dynamic formation of new vortices in a turbulent burst increases the 
vortex number immediately to A^eq- 

The evolution following a triggered injection event can be monitored by recording 
the NMR absorption height as a function of time either at the location of the 
counterflow peak, as is done in Fig. 1201 to measure the reduction in the macroscopic 
counterflow, or in the region of the Larmor edge, as is done in Fig. 1111 to exhibit a 
signal which is generated more directly by vortices. In both regions the NMR line 
shape arises from the interplay of vortex-free counterflow and of vortices via their 
orientational influence on the order parameter texture, as explained in the Appendix. 

The principal new information in Fig. [50] is related to the motion of the vortices 
after the turbulent burst. From the time interval between the trigger and the first 
response from the propagating vortex front one can measure the longitudinal velocity 
of the fastest vortices [130j. It is found to be approximately Vz ~ a^R. The subsequent 
time interval, during which the counterflow peak decays to zero, is at low temperatures 
a measure of the thickness of the vortex front: the convolution of the moving vortex 
front with certain width through a detector coil of given length. In Fig.|5n]at relatively 
high temperatures no stable front is formed. Here the time interval between the trigger 
and the moment when the counterflow peak vanishes measures the flight time of the 
slowest vortices from the AB interface to the far end of the detector coil. These 
properties of the front and the helical vortex bundle behind it will be discussed in 
Sec. |S| We note that the response in Fig. 1201 is deterministic and fully reproducible 
from one measurement to the next, in spite of the fact that it involves a short-lived 
stochastic turbulent burst. The reproducibility is attributed to the well-behaved value 
of the KH critical rotation velocity fie (T, P) plus to the fact that the propagation of 
the vortices along the rotating column controls the slow time scale in Fig. EOl which 
masks all other faster processes. 

KH injection produces the transition to turbulence at the highest onset tem- 
perature. The reason is interpreted to be the initial configuration in which the seed 
vortices end up on the B-phase side of the interface: As sketched in Fig. [21 they form 
a tightly packed bundle of many roughly parallel vortex loops, as estimated in Fig. 1131 
Apparently, in the applied B-phase flow they immediately start interacting turbulently, 
when the amplitudes of Kelvin wave excitations are not over-damped. This leads 
instantaneously to a turbulent burst and gives the highest onset temperature Ton 
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Figure 20. NMR measurement of KH instability in the turbulent temperature 
regime. (Top right panel) NMR absorption line shapes in ^He-B. The initial 
vortex-free state displays a large maximum which is shifted from the Larmor 
position by a temperature dependent amount. The height of this counterflow peak 
is roughly proportional to the counterflow velocity. The final equilibrium vortex 
state has a flat distribution of absorption. At fixed temperature the integrated 
absorption under the two line shapes is equal. (Main panel) Counterflow peak 
heights measured as a function of time in the BAB configuration, when the KH 
instability is triggered in the setup of Fig. 171 (with /j, = 8.0 A). Seed loop injection 
at the AB interface is followed by a turbulent burst, rapid polarization, and the 
forming of a propagating vortex state consisting of a front followed by a twisted 
bundle. The counterflow peak height starts decreasing when the first fastest 
vortices in the front reach the closer end of a detector coil. The peak height 
vanishes when the last vortices in the front pass through the far end of a detector 
coil. In this example at 0.55 Tc no stable front exists: the fastest vortices move 
at ~ 2mm/s while the slowest vortices, on an average, are only half as fast 
by the time they have passed through the coil. Since the measuring setup in 
Fig.|7]is symmetric, the top and bottom detector coils display identical responses. 
(Bottom panel) Externally controlled trigger for KH injection. The rotation drive 
is suddenly increased by AQ 0.05rad/s so that f2 jumps above the critical KH 
velocity of Qc = 1.48rad/s. This is signalled as an incremental increase in the 
counterflow peak heights of both signals in the main panel. 
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Figure 21. Sketch of the initial vortex loop configuration in KH injection. A 
section of the cylindrical two-phase sample is shown. Roughly the equilibrium 
number of double-quantum vortices exists on the A phase side. At the AB phase 
boundary the A-phase vortices curve into a surface layer which coats the interface 
on the A-phase side. In this layer the vortices run radially out to the cylinder wall 
IStil . The B phase is vortex free, except for the seed loops which have escaped 
across the AB interface in a KH instability event. Initially they are contained 
in a tight bundle which is oriented radially in the cylinder, starting from R^g 
at the AB interface and ending at the cylindrical wall. The bundle has roughly 
a diameter ~ 0.1 - 0.2mm (corresponding to A/2 in Fig. ll'M center), a length 
^ 0.4mm (corresponding to R — R^g in Fig. Illi . and contains on an average 10 
vortex loops with a spacing ~ 15 - 25 /im (Fig. [T^ right). In the turbulent 
temperature regime the vortices within the bundle are rapidly destabilized by 
Kelvin wave excitations, owing to the applied flow v = Vn — Vs = f2xr and the 
interactions between the loops. 



which, at least at higher flow velocities, turns out to be independent of the flow 
velocity. In contrast such transitions to turbulence, which have to be preceded by the 
single-vortex instability, depend on the flow velocity and occur at lower temperature. 

Neutron absorption: The nuclear capture reaction of a "^He nucleus with a thermal 
neutron provides an externally controllable mechanism for vortex line injection in 
vortex- free flow of ■^He-B jl31l 1^. A thermal neutron incident on liquid "^He has a 
short mean absorption length of only ~ 0.1 mm before it suffers the capture reaction 
n+'^Hc p+'^H+764keV. The reaction energy is released in the form of kinetic energy 
of the two reaction products. It is dissipated by them in ionizations and recoil, such 
that roughly a volume of the fluid within a radius < 50 /im from the reaction site is 
locally heated above Tc. Within microseconds this "neutron bubble" cools back to 
the ambient bath temperature, but a random vortex tangle is left behind TTT . In 
the absence of an externally applied counterflow the loops in the tangle shrink and 
disappear, but in vortex-free rotation the largest loops with proper orientation and 
polarization are extracted from the tangle and expand into the bulk fluid, where they 
initially appear as separated rings. 

The extraction of rings from the tangle is governed by the magnitude of the 
applied vortex-free flow velocity v according to the well-known formula for the 
equilibrium state of a vortex ring: A ring of radius r is in stable state at the flow 
velocity 

v{r) = In ^ . (57) 
^ ^ Awr ^ ^ ^ 

A ring with larger radius than r{v) will expand in the flow while a smaller will contract. 
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Therefore a minimum threshold velocity Vcn exists at which the first vortex ring can 
be extracted from the tangle. This velocity corresponds to the maximum possible ring 
size, which has the radius i?b of the neutron bubble: r{vcn) ~ ^^b- At larger flow 
velocities smaller rings can be pulled from the tangle. Simultaneously the number of 
such smaller rings can be larger than one, since several smaller rings can fit within 
the neutron bubble. Again there exists a minimum threshold velocity Vcni which is 
required in order to extract i rings of equal size from the bubble. Their number i is 
obtained from a volume argument, i.e. according to how many spheres of radius r{v) 
can fit inside the neutron bubble without overlap: i ^ {Rh/r{v))^ . Using Eq. (|57|l 
we then obtain Vcni ~ i^^^ Vcn- A better experimental and theoretical justification of 
these features can be found in Ref. [1331 IS). 

A neutron absorption event can be used to inject vortex rings in the rotating 
flow. The number and size of the rings depends on the local velocity of flow. Since 
the mean absorption length is only ~ 100 /im and the flow velocity increases with 
radius as v = fir, all rings will initially be located in the vicinity of the outer wall at 
a well-defined flow velocity v — QR. In this situation neutron absorption becomes an 
externally controllable injection method. However, because of the random tangle from 
which the injected rings originate, vortex formation from a neutron absorption event 
is by nature a stochastic process. This means that at a flow velocity flR ~ Vcni the 
number of rings obtained from a given neutron absorption event can be anything from 
zero up to the maximum limit i. At high temperatures T > Ton, each extracted ring 
evolves independently to a rectilinear vortex line. These can be individually counted 
in similar fashion as in Fig. 1111 

In the turbulent temperature regime either no vortices are obtained or, if even 
one ring is extracted from the tangle, the event leads to dynamic vortex formation and 
a turbulent burst. By choosing the value of the applied flow velocity QR in relation 
to the threshold velocity Wcni (which is measured at T > Ton), one can study the onset 
Ton as a function of the number of rings i which are extracted from a single neutron 
absorption event. Of particular interest is to check whether with flR > Vcn a single 
vortex ring can lead to a turbulent burst. In practice, because of the small size of 
the neutron bubble, i?b ^ 50 /zm, the applied flow velocity for vortex injection via a 
neutron absorption event has to be relatively high: ficn = Vcn/R ^ 1.4rad/s. Such 
measurements confirm that at sufficiently low temperatures even a single extracted 
ring leads to a turbulent burst. 

Other injection techniques: There exist also other types of "injection", in which a 
small number of curved seed vortices can be introduced in externally applied rotating 
counterflow. Vortex multiplication can, for example, be initiated by a dynamic 
remnant vortex which is left behind from a preceding experimental run and which 
has not yet had sufficient time to annihilate at the container wall A further 

source of curved seed vortex loops are those vortices of the rotating equilibrium vortex 
state which connect to the cylindrical side wall owing to the misalignment between 
the rotation axis and the symmetry axis of the cylinder |134[ I135| . The common 
feature of all such experiments with different injection mechanisms and varying initial 
conditions is that vortex multiplication and the onset of turbulence become possible 
at the temperature corresponding to g 1, although the measured onset Ton may be 
lower and depends on the initial conditions, i.e. on the applied counterflow velocity 
and the configuration and density of the initial seed vortices. 

In general, it is found both from experiments and numerical simulations that, at 
sufficiently low q, injection of even a single vortex loop into externally applied flow 
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Figure 22. Principle of the boundary induced vortex instability. A vortex ring 
of radius 0.5 mm is placed in flow with constant normal fluid velocity and its 
evolution in time is followed. Initially the center of the ring is 0.5 mm above a 
plane solid boundary, which is the xy plane at z = 0, and the plane of the ring 
is slightly tilted from the y = plane. The ring first expands and reorients itself 
in the applied flow and then collides with the boundary. The collision with the 
boundary creates two reconnection kinks on the vortex. The kink on the right has 
the correct orientation and helicity and starts to expand in the applied flow. The 
kink grows to a loop which in turn reconnects at the boundary creating thereby 
a new independent vortex. The calculation uses the hydrodynamic parameters of 
■^He-B at 0.40 Tc (g = 0.21), with the normal fluid velocity Dn = 3mm/s oriented 
along the positive y-direction. 



always leads to turbulence. Based on simulations in the single-vortex regime, where 
inter-vortex interactions can be neglected, one can conclude that the most important 
mechanism for generating a second independent vortex from the seed loop involves 
a reconnection of the seed loop at the solid boundary. The mechanism is illustrated 
in Fig. 1221 which shows a vortex ring drifting in applied flow, colliding with a plane 
boundary, and the ensuing formation of an expanding new loop adjacent to one of the 
reconnection kinks. Reconnection with the boundary induces a spectrum of helical 
Kelvin wave excitations on the vortex |13b[ll37| . At finite temperature in the absence 
of external flow these are damped, but in applied counterflow a Kelvin wave with 
correct helicity and proper orientation can start to expand jllMI . This expanding 
loop can then reconnect again with the boundary, leading to the generation of further 
loops. 

5. Propagating vortex front and twisted vortex state 
5.1. Introduction 

In the rotating container, the multiplication of the injected vortices and the ensuing 
turbulence are transient phenomena which ultimately lead to the establishment of 
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the stable equilibrium vortex state, with the equilibrium number of rectilinear vortex 
lines. In this state the superfluid component imitates solid-body rotation. Thus 
the superfluid and normal components are finally locked to corotation. In the long 
cylinder, the transient phase acquires new features, which consist of the spirally 
winding propagation of the vorticity into the metastable Landau state and the 
subsequent relaxation of twisted vortices to rectilinear lines. 

Consider a long rotating superfluid column initially in a metastable state of 
vortex-free flow with high kinetic energy: the superfluid component at rest, and 
the normal component in solid-body rotation. If nearly the equilibrium number of 
vortices is suddenly created locally in some part of the sample, how does the vorticity 
spread over the rest of the sample to reach stable equilibrium? In other words, how 
is the initially stationary superfluid component dragged into corotation? To shed 
light on this new hydrodynamic problem, we discuss here an experiment where vortex 
propagation in superfluid ^He-B is studied by monitoring the NMR signal as a function 
of time, at a location which is far away from the injection site. 

As discussed in the preceding chapter, the dynamics of quantized vortex lines 
in '^He-B is strongly influenced by the strength of the dissipative mutual- friction 
force, imposed by the normal component. This force drives the longitudinal motion 
of vortices along the rotating superfluid column, and is responsible for dissipating 
the excess kinetic energy of the initial Landau state. The reduction in this frictional 
damping with decreasing temperature is responsible for the transition to turbulence 
and is also expected to influence further the propagation on approaching the lowest 
temperatures below 0.3 Tc- 

Initially it came as a surprise that both numerical simulations and measurements 
did not bear any evidence of an expanding turbulent vortex tangle in the rotating 
sample. In fact, even under the conditions of fairly low damping (at temperatures of 
~ 0.3 Tc), the injected vorticity does not spread in the form of an incoherent tangle. 
Instead, vortices propagate along the sample in a time-invariant configuration with 
a narrow vortex front. Behind the front the vortices are left in a helically twisted 
configuration |138j . This novel dynamic vortex state arises as a consequence of the 
spiral motion of the vortex front which, in addition to the axial propagation, also has 
an azimuthal velocity component with respect to the frame of the container. The 
latter is derived from the reactive component of mutual friction. Experimentally the 
twisted vortex state can be identified through its associated superfluid velocity field 
which has a component along the rotation axis and leaves a clear fingerprint on the 
NMR signal. Both of these new features, the vortex front and the twisted vortex state, 
are also reproduced in numerical simulations. 

5.2. NMR response from propagating vortices 

We begin by examining NMR signals from vortex-injection measurements with the 
Kelvin-Helmholtz technique. As described in Sec. 14.61 vortex injection into originally 
vortex- free flow at low enough temperatures leads to a turbulent burst, where the 
number of vortex lines momentarily increases locally to iVoq in some cross section of 
the sample. The newly created vortices begin to propagate along the rotating column, 
gradually replacing the metastable vortex-free state of large rotational counterflow 
with vorticity. Figure |2S1 shows the measured NMR signals as a function of time, as 
recorded with spectrometers located at a distance of ~ 4 cm from the injection site at 
the AB phase boundary. 
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Figure 23. Time sequence of NMR absorption signals. Similar to Fig. 1201 the 
signals show the time evolution when vortex loops are injected into vortex-free 
flow from the AB interface with the setup of Fig. |3 in BAB configuration. The 
two traces represent the counterflow peak height (CF) and a maximum close to 
the Larmor edge of the NMR signal (Larmor, (•)). The corresponding frequency 
shifts are denoted with arrows in the NMR absorption spectra in the insert. The 
Larmor signal is recorded with the bottom spectrometer and the CF signal with 
the top spectrometer. They are not exactly comparable in amplitudes, but are 
normalized such that the integral over the absorption spectrum is unity. The 
spectra in the insert are measured with the bottom spectrometer as a function of 
the NMR field which has been converted to frequency shift. The Larmor region is 
recorded by sweeping the NMR field linearly back and forth in a narrow interval 
and by plotting the peak height. The location of this maximum is not exactly 
constant during the time when the Larmor absorption is time dependent. The 
CF signal is recorded at fixed NMR field. The initial state is vortex free rotation 
at 1.45rad/s. Vortex injection is triggered by a small increase of 0.05rad/s which 
is seen as a small step increase in the CF signal. At the moment of injection 
several vortex loops cross the phase boundaries, undergo a burst of turbulence, 
and start the propagation towards the detectors at both ends of the sample. No 
change in the signals is seen until the first vortices reach the closest end of the 
detector coil and the CF signal starts decreasing owing to the removal of the 
azimuthal counterfiow. Following this, the Larmor signal rises to a maximum 
hmax which provides a measure of the axial counterfiow velocity in the twisted 
vortex state. Subsequently the Larmor signal decays exponentially in amplitude 
which corresponds to the relaxation of the twist when the vortices are connected 
at one end to a solid end plate of the sample cylinder. 
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The NMR signals in Fig. ESI represent traces of the absorption as a function of 
time at two different locations of the NMR spectrum, corresponding to two distinct 
maxima in the NMR line shape. The relative amount of absorption concentrated 
at these peaks is a useful measure of the azimuthal counterflow and, consequently, 
of local superfluid vorticity: if the counterflow is small, an absorption peak forms 
near the Larmor frequency, whereas in situations with large counterflow most of the 
absorption collects at another peak at slightly higher frequency, which corresponds to 
the so-called counterflow peak (see insert of Fig. I23|l . In the main panel the heights 
of the Larmor and counterflow peaks (marked as CF) are recorded as a function of 
time. Actually, in the precise setup shown in Fig.|7|the two peaks are being monitored 
with two independent spectrometers located near the opposite ends of the cylindrical 
sample. In KH injection these two NMR traces refer to different, disconnected sets 
of propagating vortices, but since the sample arrangement is exactly symmetrical in 
the two halves, the two signals can be viewed as representing NMR responses on a 
common time scale. 

At the start of the experiment, the B-phase volume contains no vortices and is in 
a state of high counterflow at = 1.45 rad/s: a large CF signal is then visible, whereas 
in the Larmor region there is very little absorption. At time t = 0, the angular velocity 
is suddenly increased to 1.50 rad/s (seen as a small increase in the CF peak height), 
above the KH critical velocity which has been adjusted between 1.45 and 1.50rad/s 
with the barrier magnetic field. At this moment, vortices are injected into the B phase 
sections, undergo the turbulent burst, and begin to propagate towards the NMR coils. 
After a temperature-dependent flight time (approximately 20 s in the conditions of 
Fig. I23|l the vortices arrive at the two coils, which is signalled by a rapid decrease in 
the measured azimuthal counterflow and eventually by a complete disappearance of 
the CF signal. Simultaneously, the peak height of the Larmor signal increases, reaches 
a maximum, and then slowly relaxes to the value corresponding to the equilibrium 
NMR spectrum, with roughly the equilibrium number of rectilinear vortex lines. 

Two interesting conclusions can be immediately drawn from the responses in 
Fig. First, the fact that after the vorticity arrives at the NMR coil the 

counterflow signal drops rapidly from its initial maximum value to zero indicates 
that the propagating vortices form a front, which separates the vortex- free region 
from the region occupied by vortices. Similar behavior is observed also in numerical 
calculations, which will be discussed below in Sec. 15.61 

The second peculiar feature has to do with the response of the Larmor signal, 
its nonmonotonic time dependence. The transient overshoot of the Larmor peak 
appreciably above the equilibrium-state value is especially noteworthy. Such a 
response cannot be obtained with any configuration which would incorporate only 
rectilinear lines. Instead, this is a fingerprint of the helical twist acquired by the 
propagating vortices. In the following, we discuss this aspect in more detail and 
demonstrate how both numerically calculated NMR spectra for the twisted vortex 
state and detailed simulations of vortex motion support this interpretation. 

5.3. Helically twisted vortex state 

To understand the NMR features discussed above it is instructive to begin by studying 
the motion of a single vortex loop in rotating flow. Consider the instantaneous velocity 
of the vortex-line element s at the point where the line connects to the cylindrical side 
wall in Fig. 1241 In the rotating frame, the velocity of the normal component vanishes. 
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Figure 24. The motion of a single vortex loop in a rotating cylinder. If viewed 
in the frame rotating with the container, the end of the loop that connects to the 
side wall moves with a velocity that has both an azimuthal as well as an axial 
component. Therefore, the plane of the loop rotates while the vortex expands in 
the rotating flow. The rest of the vortex is aligned along the central axis which 
is the equilibrium position of the vortex. 



while the local superfluid velocity is Vg = —ilR 4> (here we ignore the small self- 
induced contribution from the curvature of the loop). From Eq. one then finds 
the velocity of the element s as 

^i^^ -{I- a')Q.R ^- aVLRi, (58) 

i.e. the vortex end has both azimuthal and vertical velocity components, as shown 
in Fig. 1^ From this, we immediately find that the expansion time, or the "time of 
flight" which is required for the single curved vortex to advance an axial distance d 
along the rotating column, is d/{aVlR). Combined with this, the azimuthal velocity 
gives rise to a spiralling motion where the plane of the loop is constantly turning 
during its expansion. 

Consider next the situation with a large number of loops (close to iVcq) in a 
configuration similar to that in Fig. |31 one end connecting to the top end plate, and 
the other to the side wall. The local superfluid velocity at any element on one of these 
loops is then modified by the contribution induced by the other vortices. In particular, 
the vortex segments at the top end plate would be expected to have Vg ~ Vn, and be in 
solid-body rotation with the container. Hence the azimuthal counterflow experienced 
by the vortex elements changes from zero to ~ fir over the volume occupied by 
vortices. On the other hand, both numerical simulations and experimental data (see 
below for more details) indicate that the flight time for a cluster of vortices remains 
almost unchanged from the single-vortex value. The combined effect of these two 
motions is to drive the vortex cluster into a helically twisted configuration, as shown 
schematically in Fig. 123 Estimating (1 — a')J7 as the angular velocity of azimuthal 
motion, and taking aVLR for the axial velocity, we arrive at the approximate value 
k ~ (qR)"^ for the wave number of the helical twist. 

A more extended discussion of the formation and properties of the twisted vortex 
state can be found in Ref. jl88| . Most importantly, however, the helical structure 
implies superflow in the axial direction which, as we now proceed to discuss, explains 
the peculiar transient overshoot in the NMR absorption near the Larmor frequency, 
as observed in the measurement of Fig. 1231 
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Figure 25. The vortex front expands along the rotating column into the region 
with vortex free counterflow with the longitudinal velocity vl,z ~ af!R. In the 
NMR signal the propagating vortex front is manifested by the rapid disappearance 
of the counterflow signal and by the increase in absorption in the Larmor region 
of the spectrum (Fis. 1231 . Behind the front the number of vortex lines is close to 
that in equilibrium and also superfluid flow is much closer to the ultimate state of 
solid body rotation. The vortices in the front move also in the azimuthal direction 
with the velocity vl,cI> ~ ct'flR and thus generate a helically twisted vortex bundle. 
The twist gives rise to superflow in the z-direction with the velocity Vb^z such that 
at the center of the cylinder there is upward flow and close to the wall downward 
flow (the case in this flgure where the front propagates down). The helical twist 
with its z-directional flow is seen as an overshoot of the NMR absorption signal in 
the Larmor region of the spectrum, while the relaxation of the twist explains the 
exponential decay of the Larmor overshoot ("Fig. I2'.ii . The velocities in this figure 
are marked as viewed from the inertial laboratory frame and the flow velocities 
according to the uniform twist model. 



5.4- Superflow field of twisted state: model of uniform twist 

We will now construct a model for the superfluid-velocity profile in the twisted 
vortex state. This is achieved by noting that the equation of motion for superfluid 
hydrodynamics, Eq. H16() . allows stationary (dvs/dt = 0) solutions which represent 
twisted vorticity. The simplest one of these is translationally invariant along the axial 
(z) direction and of the form 



This situation corresponds to one where the wave vector k of the twist is independent 
of the radial coordinate r, or 



Here = u^cj) + is the local coarse-grained superfluid vorticity. The condition 
H60|l simply means that in the case of uniform twist the azimuthal tilt of the vortices 



Vs = v^{r) 4> + Vz{r) z. 



(59) 




(60) 
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v.ir) ^ " , )' ■ (62) 



has to increase with radial distance from the cyhnder axis. In addition, we impose 
a further condition on Vg by requiring that the mutual-fr'iction-induced force on the 
vortices vanishes, i.e. 

LJs X (Vs - Vn) = 0, (61) 

meaning that the vortices are ahgned paraUel to the counterfLow in a "force-free" 
configuration. An equivalent condition would be to require that radial vortex motion 
is absent: vi^^r = from Eq. Under these requirements, the superfluid velocity 
field can be solved as 

fir + {kr)vo 

= l + (fc,)2 - 

Vq — (fcr)f2r 
1 + (fcr)2 

The parameter vq is fixed to Ug — (ri/fc){(fci?)^/ ln[l + (fci?)^] — 1} by requiring that 
the net flow through any cross section of the cylinder vanishes. 

The NMR response from a twisted vortex state with a flow profile like that in 
Eq. 1)62(1 can be roughly understood in terms of the following simple arguments. The 
local resonance frequency depends on the orientation of the B phase order parameter 
which, in turn, is affected by the magnitude - and orientation - of the local counterflow 
velocity. This can be described in terms of a free-energy contribution which has the 
form —[1b • (vn — Vg)]^, where 1b = [— z + 5(z • n)n -I- vT5(z x n)]/4 is the orbital 
anisotropy axis (in the presence of an axial magnetic field) corresponding to the order- 
parameter vector n. Therefore, counterfiow along the cylinder axis should have a 
tendency to favour the orientation 1b || n || z. This in turn corresponds to additional 
NMR absorption at the Larmor frequency, as observed in Fig. As shown below 
in Sec. 15.51 these arguments are further supported by more detailed comparison to 
order-parameter textures which are associated with flow profiles of the form given by 
Eq. (j^. 



5.5. Experimental results on twisted vortex state 

Next we note some further experimental details on the propagating vortices and the 
twisted vortex state. Fig. shows measurements on the flight time of the vortex 
front as a function of temperature. In these experiments, the Kelvin-Hclmholtz 
instability is triggered and the time of flight ^flight for the vortex lines to reach the 
closer end of the pick-up coil (a known distance d away) is recorded. The data has been 
presented in terms of the mutual-friction parameter a, as extracted from the relation 
a = d/ {tfiight^R), which is based on the assumption that the vortices propagate with 
axial velocity afii?. The degree of validity of this assumption can be assessed by 
comparing to the previously published data of a(T) from Ref. 72 . 

The division of vortex propagation into regular and turbulent regimes as a 
function of temperature means that different experimental procedures had to be 
followed, to achieve as similar initial conditions for the measurements as possible. 
In the low-temperature (low damping) regime, T < 0.6 Tc, a KH vortex injection 
always initiates turbulent vortex formation, where nearly the equilibrium number of 
vortices {N « TVeq) is created and then propagates along the column. However, in the 
regular regime at temperatures above 0.6 Tc the vortex number is conserved and only 
the few initially injected loops expand in the column. In this regime a much larger 
number of vortices can be injected using the Kelvin-Helmholtz method differently: the 
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Figure 26. The dissipative mutual friction coefficient a, plotted on a logarithmic 
scale as a function of temperature. It is determined from the axial flight 
time, assuming «l z = ceflR, regardless of the number of vortex lines in the 
front. Flight time data have been used from NMR measurements to extract 
a = d/ (tfiigi^^QR) where d is the distance which the vortex lines have traveled. 
These are marked with large size filled symbols, which correspond to different 
methods of determining the flight time (for details see TW). The data marked 
(a) is measured using the method explained in Figs. l2(11 and l23l The small-size 
symbols (•) represent the data from Ref. \72\ for comparison. The agreement 
in a values between the different sets of data is well within the accuracy of the 
individual measurements, which can be argued to justify the model for the flight 
time of the first vortices in the vortex front at T > 0.45 Tc- 



sample is rotated at high > 2 rad/s (weh above the KH critical velocity) with only 
B phase in the column. Next the barrier field is increased by an incremental amount 
so that suddenly a narrow sliver of A phase is created in the middle of the column. 
When the A phase forms, the AB interfaces are unstable from the start, until a large 
number of vortices have been transferred to the surrounding B phase. 

As mentioned above, the characteristic signature of the twisted vortex state in the 
NMR spectrum is the absorption overshoot in the vicinity of the Larmor frequency. 
A quantitative measure is the maximum peak height of the overshoot during its time- 
dependent evolution, which is denoted with the symbol ft-max in Fig.|231 Fig. 1271 shows 
the dependence of the maximum peak height on temperature. Here the peak height 
is normalized to the corresponding peak height ho of the NMR spectrum from the 
nonrotating state at Q = 0. The solid line in the figure represents the corresponding 
quantity obtained from numerically calculated NMR line shapes. In these calculations, 
the equilibrium order-parameter distributions were determined by minimizing the 
textural free energy in the presence of superflow according to the uniform-twist model, 
Eq. (|62(l . The wave vector of the twist was chosen as fc = 1/ {qR) (see discussion in 
Sec. 15.3(1 . Towards low temperatures the temperature dependence of the Larmor 
overshoot is seen to reflect the exponential dependence of on 1/T: the twist 
becomes tighter with decreasing temperature and leads to an increased superfluid 
velocity in the axial direction, and hence to more excess absorption near the Larmor 
frequency. The dependence of the overshoot on the angular velocity is presented on 
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Figure 27. (Left) The amplitude of the maximum Larmor overshoot /imax : 
divided by the maximum hg of the corresponding absorption spectrum in the 
nonrotating state £7 = 0, plotted as a function of temperature. The measurements 
have been performed at a pressure of P = 29.0 bar. Below 0.60 Tc the rotation 
velocity f2 is as shown for the KH instability in Fig. 1141 while above 0.7 Tc a 
higher rotation velocity of Q = 2.5 rad/s is used. The solid line shows the 
corresponding quantity from numerically calculated NMR absorption spectra, 
using the uniform-twist model with kR = l/q. (Right) The amplitude of the 
maximum overshoot divided by the maximum of the f2 = spectrum as a function 
of fi. These measurements have been performed in the temperature interval 0.50 
- 0.55 Tc. The solid line is the theoretical prediction of the uniform-twist model 
with kR= l/q. 



the right m Fig. |23 

The theoretical curves in Fig. have been calculated without adjustable 
parameters. In view of this, the agreement between experiment and theory can be 
considered as remarkably good at higher temperatures where the model of uniform 
twist could be expected to apply. Thus the calculations of the NMR spectra 
corroborate the explanation of the NMR measurements in terms of a traveling vortex 
front and a trailing twisted vortex state. Since the '^He-B NMR absorption spectrum 
does not uniquely define the order parameter texture (see Appendix), the agreement 
in Fig. 1^ might not entirely exclude other possibilities. Nevertheless, it appears safe 
to conclude that a new dynamical vortex state has been identified. In the following 
section we describe numerical simulations of the propagation of vortices in the rotating 
column. They provide additional justification for both the vortex front and the twisted 
state. 

5.6. Propagating vortex state in simulations 

Our numerical simulations are performed using the vortex filament model |139l I114| 
where vortices are included individually as topologically stable line objects. The 
local superfluid velocity at a particular vortex segment s is calculated from a Biot- 
Savart integral by summing the flow contributions from all other vortex elements. A 
discretized Laplace equation is solved to take into account boundary conditions on 
vortices and their image fields in the geometry of a finite cylinder (for details see 
Ref. [HOI). 

Consider first the case of a single vortex expanding in rotation fFig. l24|l . Ignoring 
the vortex curvature, the radially oriented endpoint on the cylinder wall propagates 
along the cylinder with the velocity aClR and rotates with angular velocity a'Q (in the 
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Figure 28. Vortex configuration (left), longitudinal counterflow velocity V2,(r) 
(center), and azimuthal counterflow velocity v^{r) (right) at t = 60s after the 
start of the propagation along a cylinder of radius R = 1.5 mm and length L = 
40mm at 0.40Tc and Q = l.Orad/s. Only a 30mm long section of the lower 
end of the cylinder is shown. For clarity the axial and radial dimensions have 
diff'erent scales in the left panel. The velocities in the center and right panels 
are shown at radii r = nR/6, with integer n. The spikes on the r = traces 
originate from the noise in the position of the center-most vortex with respect to 
the axis of the cylinder. The boundary condition on the vortices at the bottom 
end plate of the cylinder at z = causes the twist and Vz{r) to vanish there. A 
stable time-invariant vortex front travels along the column and has here reached 
a height z ^ 12 - 14 i?. This is evident from the fact that the radial velocity 
profiles Vz{r) and v^(r) do not change as a function of time t. 



laboratory frame) . The self- induced superfluid velocity of the curved vortex (together 
with the boundary contribution due to the image field) is mainly along Vn and reduces 
the counterflow velocity. This does not affect the longitudinal propagation velocity 
substantially if the rotation velocity is large, but at low H. the difference is easily 
observable, and a somewhat slower propagation speed than aQR is realized. Similarly, 
the self-induced velocity tends to speed up the azimuthal rotation of the vortex. This 
is especially visible at low temperatures where a' is small and the self-induced velocity 
gives the main contribution to the rotation. A smaller cylinder radius enhances these 
effects further. 

The same principles are at work when many vortices expand along the column. 
Such calculations are started by placing N^q curved seed vortices in the cylinder. The 
usual starting configuration is that with one end of the vortex on the bottom end 
plate of the cylinder and the other end bent to the cylindrical side wall. The locations 
of the vortex ends on the bottom plate are the equilibrium positions of rectilinear 
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vortex lines in solid-body rotation. The choice of these locations and the way how the 
vortices are bent to the cylindrical wall prove to have a minor effect on the outcome 
from the simulations. This is because during the first few seconds the vortices move 
rapidly around owing to their strong curvature, which upsets, for instance, the ordered 
initial configuration on the bottom end plate. 

In the measurements A'oq is of order 10"^. For numerical calculations such a 
large number of vortices is too time-consuming, since the spatial resolution has to be 
kept below the average inter-vortex spacing (~ {2Q/ k)~'^^'^ ~ 0.2mm). In practice 
this limits the calculations to low rotation velocities or, as was done here, to smaller 
cylinder radii than in the experiments. With a large number of vortices there are 
more reconnections between vortices, which partly disrupts and straightens the twisted 
state. Figs. 1281 and 1291 illustrate the characteristics of the propagating front and the 
twist which the former generates at 0.40 Tc and 0.60 Tc temperatures. At the lower 
temperature of 0.40 Tc a sharp and stable vortex front is formed, but at the higher 
temperature of 0.60 Tc this is not the case. Thus the stability of the front is strongly 
temperature dependent, although in both cases the trailing vortex bundle behind the 
front is clearly twisted. 

The stability of the front is governed by the axial counterflow velocity and thus 
by the pitch of the twist immediately behind the front. At high temperatures the 
axial counterflow is weak. Here the vortices which happen to fall behind the front 
feel a reduced azimuthal counterflow Vn — Ws,0 fFig. I29l (right)) and, as a result, their 
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Figure 30. (Left) Vortex-front velocity in simulations with a cylinder of radius 
R = 1.5mm and length L = 40mm rotating at f2 = l.Orad/s (expressed in the 
rotating frame). The axial velocity f front z (denoted with o) was obtained from 
the location z of the front at different points in time. The front was defined to 
be at the location z where the azimuthal counterflow at r = i? had decreased 
2% from its maximum value. The azimuthal velocity ffront,^ (denoted with x) 
is the average velocity of the vortex endpoint at r = R when the front is at 
z Si L/2. In the averaging only the vortices within 5 mm from the front were 
taken into account. Note that the enhancement of the axial velocity above the 
value Dfi-ont.z ~ aflR at temperatures below 0.45 Tc is caused by the longitudinal 
velocity contribution from the twist. (Right) Longitudinal width or thickness of 
the front in the same simulation calculations. At low temperatures below 0.45 Tc 
the tight twist generates a large longitudinal velocity contribution which stabilizes 
a narrow front. At high temperatures vortices lag more and more behind the front 
and its width slowly increases during the propagation along the cylinder. In this 
regime the propagating structure does not qualify as a time invariant traveling 
solution. 



longitudinal propagation slows down. This is seen in simulations at high temperatures 
where the longitudinal spread of the vortex ends on the cylindrical wall increases 
continuously with time. At lower temperatures the axial flow is larger, reaches a 
clear maximum just behind the front (Figs. l^ ( center)), and provides a compensating 
longitudinal force on the vortices which fall behind. This is seen by considering the 
motion of a radially oriented vortex segment ending on the cylindrical wall. The 
longitudinal velocity of this vortex element is given by wl.z = o:{vYi — Vs.^) + (^l — a')vs^z- 
On moving towards low temperatures, where both a and a' approach zero, but are 
similar in magnitude, the increasing axial velocity Vs^z behind the front is sufficient to 
compensate for a reduced azimuthal counterflow, if a vortex falls behind the front. 

Above 0.45 Tc the width of the front increases with time. This is seen in the right 
panel of Fig. 1201 which illustrates the front thickness as a function of its axial location. 
A time-invariant propagating front solution is only obtained at temperatures below 
about 0.45 Tc. Thus the definition of the front at temperatures above 0.45 Tc becomes 
arbitrary. However, to characterize the propagation we continue assigning widths to 
fronts up to 0.7 Tc in FiglJOl We locate the position of the front to where the azimuthal 
counterflow at r = i? has decreased 2% from its maximum value flR. The rear of the 
front is chosen to be located at the position where the azimuthal counterflow at r 
= R falls by one half of the difference between its maximum and minimum values 
(the latter is nonzero because of the vortex-free region outside the equilibrium vortex 
bundle). We then define the front width as the difference between these two positions. 
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In the left panel of Fig|2II|the longitudinal and azimuthal propagation velocities of 
the vortex front are shown at different temperatures above 0.30 Tc- Above 0.45 Tc the 
longitudinal propagation velocity is approximately given by aflR but on cooling below 
0.45 Tc the speed starts to deviate from this value. The rapid rise of the normalized 
propagation velocity t)front,z/ (aOi?) in this regime is a clear indication from the strong 
pushing action which arises from the axial counterflow, generated by the tightening 
twist. This enhancement in the longitudinal velocity agrees with recent measurements 
at temperatures below 0.40 Tc. 

At high temperatures above 0.7 Tc the twist is barely visible, but it grows rapidly 
tighter with decreasing temperature. One of the evident difficulties in the present 
interpretations is the stability of the twisted state on cooling to temperatures below 
0.30 Tc. The wave vector of the twist is approximately proportional to 1/g = (1— a') /a. 
Ultimately at the lowest temperatures the helical vortices may become unstable with 
respect to Kelvin wave formation, which will lead to vortex reconnections and vortex 
multiplication. The details of such breakup of the front depend on temperature, vortex 
number, and rotation velocity. With R — 3 mm, fl — Irad/s and with 23 vortices 
the twist is stable at 0.40 Tc, while at 0.30 Tc new vortices are generated owing to the 
break up in the helical structure by vortex reconnections, if the original number of 
vortices is much less than in equilibrium. From these examples it becomes evident 
that the spiralling propagation of vortices in a rotating column cannot be extrapolated 
far below 0.30 Tc before other phenomena can be expected to intervene. Therefore 
it becomes more and more pressing to understand the correct behavior for vortex 
propagation on approaching the zero temperature limit. A second related problem is 
the fact that in numerical calculations vortices connecting to the cylindrical side wall 
of the rotating column tend to be more stable than in experiments. The origin for 
this difference has not yet been explained. 

6. Concluding remarks 

We have discussed new developments in superfluid hydrodynamics, which have been 
identified from measurements on "^He superfluids. These measurements pertain to 
situations where the velocity of the externally applied counterflow dominates over 
other contributions, which might arise from other existing vortices. The applied flow 
is generated with uniform rotation. This is chosen for experimental convenience, it 
does not dictate the existence of the observed phenomena, except for the twisted 
vortex state which is a true peculiarity of rotation. Partly the new phenomena appear 
owing to the multicomponent order parameter structure of the '^He superfluids. Such 
effects cannot be reproduced with the traditional ^He-II superfluid, where only the 
U{1) symmetry is spontaneously broken. Other phenomena were revealed because of 
the different measuring techniques which are employed in superfluid '^He research, in 
this case uniform rotation of a long cylindrical sample with multiple detectors. The 
common link, which connects the different observations, is the possibility to stabilize in 
the long sample one or two AB interfaces in a two-phase arrangement with vortex-free 
counterflow in '^He-B. 

The starting point is the superfluid Kelvin-Helmholtz instability of the AB 
interface, which becomes possible owing to the difference of three orders in magnitude 
in the vortex core radii of '^He-B and ■^He-A. In the former case quantized vortices 
have more traditional structure while in the latter case continuous vorticity is generally 
formed which is composed of doubly quantized vortex-skyrmions or of a meandering 
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vortex sheet where the vortex-skyrmions have been confined as a hnear periodic chain 
within a domain- wall- hke soliton sheet. The critical velocity for the formation of 
continuous vorticity is more than an order of magnitude smaller than for singular 
^He-B vortices. This large difference makes it possible to have a shear flow state in 
the rotating cylinder, where an equilibrium vortex state of '^He-A coexists with the 
vortex-free Landau state of ^He-B across the AB interface. 

Experimentally the Kelvin-Helmholtz instability of the AB interface is an un- 
usually robust phenomenon, with a well-behaved critical velocity for the nondissipative 
shear flow state of the two superfluids. Theoretically the instability can be associated 
with two critical thresholds of different nature. The higher threshold corresponds to 
the traditional critical velocity in the relative motion of two ideal inviscid liquids, i.e. 
it follows the expression for the classical Kelvin-Helmholtz instability in the absence of 
viscosity. This threshold might be the appropriate one in the limit of zero temperature 
far from the wall of the container with no contact to a reference frame. Under the 
conditions of the present measurements at temperatures above 0.3 Tc, the contact is 
established by the equilibrium distribution of quasiparticle excitations and the AB 
interface instability occurs at lower critical velocity. This second lower threshold 
is also fundamental in the sense that it does not depend on the magnitude of the 
interaction with the environment, the only requirement is the presence of a reference 
frame imposed by the environment. In equilibrium this is the frame where the normal 
component is at rest. As a result, the lower threshold is not determined by the relative 
motion of the superfluid components, but by the relative motion of each component 
with respect to the normal component. Interestingly, both thresholds can be compared 
to the instability of quantum vacuum around the black hole. While the lower threshold 
corresponds to the instability of quantum vacuum in the ergoregion close outside the 
event horizon, the upper threshold is related to the instability within the event horizon 
owing to the black-hole singularity. 

The Kelvin-Helmholtz instability of the AB interface has become important as 
a tool to inject reproducibly many small vortex loops in tight proximity with each 
other in B-phase flow. This process was the key to the identiflcation of the transition 
to turbulence as a function of mutual friction, controlled by the intrinsic damping 
parameter 1/q = (1 — a')/a, which here is the equivalent of the Reynolds number. 
Kelvin-Helmholtz injection starts the transition at the lowest value of 1/g and is 
now believed to be the closest example to the limiting case of a velocity independent 
direct transition into turbulence where the initially injected seed vortices interact and 
immediately start a turbulent increase in vortex number. At higher temperatures 
with a lower value of 1/q no examples exist of vortex multiplication, even if large 
numbers of vortices are injected in the applied flow. In contrast, at lower temperatures 
and higher values of 1 / q the transition becomes increasingly more probable when the 
applied flow velocity is increased and the conflguration of the injected vortices is more 
favourable for interactions between the expanding seeds. Eventually at the lowest 
temperatures ~ 0.30 Tc a vortex injection of any kind will always lead to turbulence. 
Overall the most unexpected new result is the fact that in the limiting case of a direct 
transition to turbulence, like Kelvin-Helmholtz injection in Fig. the threshold 
becomes velocity independent, as can also be concluded from numerical simulations 
and theoretical models. These models now show that in the theory of homogeneous 
developed superfluid turbulence the Kolmogorov-Richardson cascade is modified by 
the influence of mutual friction. 

In ^He-B the transition to turbulence at 1/q ^ 1 is in the middle of the 
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experimental temperature range in stark contrast to "^He-II. This is a property 
expected for Fermi liquids in general. In most superfluid and superconducting systems 
with Cooper pairing it is to be expected that g(T)~^ crosses unity at T 0.5 Tc, if the 
spectrum of fermionic quasiparticle excitations is fully gapped and the system is in the 
clean limit where impurity effects can be neglected. The underlying mechanism for 
mutual friction is in this case analogous to the chiral anomaly in relativistic quantum 
field theories (4,. 

The final topic in this review concerns the evolution of the vorticity after injection 
in the long rotating column. Turbulence in the rotating column is a short burst where 
the number of vortex lines rapidly increases. These are immediately polarized by 
the rotating flow. After the burst the vortices expand towards the unstable vortex- 
free flow in the form of a sharp front which is followed by a twisted vortex bundle. 
The expansion is governed by the mutual friction controlled motion of the vortex 
lines in the front which travel in cork-screw-like manner along the cylinder. These 
observations raise a host of new questions which await further work, especially at lower 
temperatures, to elucidate the characteristics of the T limit. Can the spiralling 
propagation of vortices along the rotating column be expected to remain stable on 
approaching the zero temperature limit? Is there some new form of dissipation present 
at the lowest temperatures, as suggested by recent measurements ^H] of the decay of 
vibrating wire or grid generated turbulence at temperatures below 0.2 Tc in "^He-B? 
It is also to be expected that studies of vortex dynamics in rotating Bose-Einstein 
condensates will shed more light on these questions |141l I142j . 

Finally we wish to emphasize that, as distinct from the traditional superfluid 
example case of ^He-II, '^He superfluids incorporate interacting Fermi and Bose 
quantum fields. This is similar to the situation in relativistic quantum field theory. 
In a multi-component order parameter field the analogy becomes stronger such that 
quasiparticle excitations become the counterparts of elementary particles (electrons, 
neutrinos, quarks), while the order parameter collective modes resemble the gauge 
fields and gravity. Such analogies provide contacts to other fields of physics and a 
strong incentive to learn more about the coherent quantum systems in condensed 
matter. They also present astonishing evidence for the underlying universality of 
physical principles in seemingly different theories. 
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Appendix: Nuclear magnetic resonance in •^He-B 

Continuous wave nuclear magnetic resonance measurement (cw-NMR) can be used 
both in '^He-A and '^He-B to monitor the state of the rotating sample and to retrieve 
information about the structure and configuration of quantized vorticity. In both 
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phases the measurement entails the task of working out the order parameter texture 
from the recorded NMR absorption spectrum, albeit in different ways. The B-phase 
spectrum provides a mapping of the spatial variation of the global order parameter 
texture. This texture is affected by the presence or absence of vortices. Thus indirectly 
from the changes in the global texture one can infer information about vortices, up to 
the point that in certain situations the vortices can be detected individually, as seen in 
Figs.r^andllll This is different from the A phase where NMR absorption is sensitive to 
the local order parameter distribution, in particular in locations where its alignments 
differ from the spin-orbit equilibrium configuration. Here each type of topological 
defect leaves a characteristic satellite peak in the absorption spectrum il43j. The 
frequency shift of the satellite identifies the type of defect and the height of the peak 
the number of these defects. Thus "spectrometry" of different types of defects in the 
A-phase order parameter texture is even more straightforward and powerful than in 
the B phase. Since B-phase measurements are in the forefront in this context, in 
the following we concentrate on the main characteristics of B-phase order parameter 
textures: the influence of vortices on the texture and how this is reproduced in the 
NMR absorption spectrum jl44| . Technical requirements for single- vortex resolution 
and cw-NMR spectrometer design can be found, for instance, in Ref. jl45j . 

The bulk tensor order parameter of undisturbed '^He-B can be written in the 
form A = A e*'^R(n, 6*), where the (real- valued) quantities A(T, P) and (p are 
the magnitude and the phase of the order parameter, and R is a rotation matrix 
which can be parametrized in terms of a rotation axis fi and an angle 9 as i?^^ = 
cos 9 dfj^j + (1 — cos 9) fi^fij — sin 9 e^jkfik- While A and 4> are usually fixed by external 
constraints and 9 by the dipolar spin-orbit interaction, the orientation of the unit 
vector n is influenced by various orientational interactions, giving rise to n(r) textures. 
The appropriate equilibrium texture is the minimum of a free-energy functional, the 
main contributions to which in the presence of an external magnetic field H are the 
field anisotropy energy 




(A.l) 



the surface energy (the unit vector s denotes the surface normal) 




(A.2) 



the energy due to the counterflow velocity field v(r) 



v„(r) - Vs(r) 




(A.3) 



and the vortex contribution 




(A.4) 



where 1 is a unit vector directed along the vortex line and the integration extends over 
the volume occupied by vortices (the information on the vortex density is contained 
in Alh)- Our notation follows that of Ref. |14()j . In addition to the above energy 
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R=:2.9mm, H=20mT, Q = 2 rad/s 
Ar= 320 lines, r=0.75Tc, P=29bar 



Figure Al. Axially symmetric flare-out n texture in an infinitely long cylinder, 
when the applied field is oriented along the symmetry axis of the cylinder. The 
arrows denote the projection of fi in the plane perpendicular to the cylinder axis. 
The shaded region in the center is the vortex cluster with 320 rectilinear vortex 
lines. The texture has been calculated with the parameters given under the cross 
section through the cylinder. 



terms, spatially varying order-parameter distributions are associated with a gradient 
(or bending) energy 



where the first integral is taken over the volume and the second over the surface of 
the sample. The characteristic length scale of the B-phase textures can be obtained 
by balancing the gradient energy with the bulk magnetic energy Fdh: the resulting 
magnetic coherence length is defined as — ^65AG2/(8ai/2), and is inversely 
proportional to the magnitude of the external magnetic field. Since the typical 
values of are of the order of a millimeter, the textures are extended and in a 
usual experimental setup the finite container size and the associated boundary effects 
become important and need to be taken into account. 

Also, it should be noted that the presence of vortices in the sample volume 
modifies the texture in several different ways. The free-energy term _Flh contains both 
the effect of order-parameter suppression at the vortex cores, and the orientational 
effect of the quantized superflow fields circulating the cores. Additionally, the vortices 
modify the global counterflow velocity field v(r) entering the expression of i^nv; this 
is typically the dominating effect. 

In a long cylindrical container with the external magnetic field directed along the 
axis, the minimum-energy n texture is axially symmetric of "flare-out" form |147II14?^ . 




(A.5) 
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Figure A2. Calculated radial distributions of the azimuthal angle a(r) (left) 
and polar angle /3(r') (right) of the B-phase order parameter orientation, the unit 
vector n, in the flare-out textures of the three NMR spectra in Fig. lASl 



In cylindrical coordinates, this can be represented in the form 

ri(r) = — sin/3(r) cosa(r) f + sin /3(r) sin Q;(r) (/} + cos (3{r) z. (A. 6) 
The NMR resonance frequency v depends on the polar angle f3 as 

^(/?)«'^o + :^sin2/3, (A.7) 

where vq = 7_ff/(27r) and v^{T,P) are the Larmor frequency and the B-phase 
longitudinal resonance frequency |149[I15(1) . respectively. This approximate expression 
is valid at high external fields, when J^o, which is the usual experimental regime. 

In the local oscillator picture (excluding line broadening effects), a spatially varying 
(3{r) leads to a distribution of resonance frequencies, which then define the absorption 
spectrum |15U) 

P{v) (X j dr r 5[v ~ v{r)], (A.8) 

where v{r) = v[(3{r)]. This defines the connection between the order-parameter 
texture fi{r), found by minimizing the sum of Eqs. ljA.l(l ~ (jA.5|l . and the measured 
line shape P{y)- 

The form of the relevant free-energy functional in rotating '^He-B is complicated, 
giving rise to several different types of equilibrium textures depending on the 
magnitude and direction of the counterflow in the sample |151II152] . However, a rough 
qualitative view of the textures can be obtained relatively easily. First, consider the 
situation in axial magnetic field and zero counterflow: through Eq. IjA.ljl the field 
tends to align n || H, or sir? (3 ~ 0, in bulk. Through Eq. (|A.7|I . this corresponds 
to NMR absorption at the Larmor frequency vq. However, the surface energy at the 
cylindrical sidewall, Eq. (|A.2|I with s = — f, favors the orientation siv? (3 = 4/5 and 
shifts the NMR absorption to higher frequencies. A combination of these two orienting 
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Figure A3. Examples of NMR absorption spectra measured witli tiie bottom 
spectrometer in the setup of Fig. 13 Tlie measurements are performed with a 
highly tuned resonance tank circuit, so that the excitation frequency is kept 
constant while the external polarizing field is swept to record the line shape. The 
Larmor field is at 21 mT. This is the value around which the NMR absorption 
is centered in the normal phase at temperatures above Tc. In the B-phase the 
Larmor value becomes a sharp edge which borders the absorption towards high 
fields (or equivalently, low frequencies). The oscillations close to the Larmor 
edge are spin wave resonances which grow larger in amplitude on cooling to 
lower temperatures. The three line shapes have been measured for (1) vortex- 
free rotation (A'' = 0, Q ^ 0), (2) equilibrium vortex state {N = A^eq), and (3) 
stationary state (nonrotating, Q = 0). 



effects and of the gradient energy creates the characteristic even distribution of NMR 
absorption in the nonrotating state of Fig. IA3I 

If the sample is rotated with angular velocity in the vortex-free state, then 
azimuthal counterflow v(r) = Or (p is created. Since the flow energy Fny in this case 
is also minimized by having sin^/3 = 4/5, this leads to the formation of the so-called 
counterflow peak in the NMR spectrum (Fig. IA3(I . In contrast, in the equilibrium 
vortex state the sample is filled with rectilinear vortices oriented parallel to the axis of 
rotation, evenly distributed with an areal density = 2n/K, such that the counterflow 
vanishes on an average. In this case, the NMR line shape is similar to that of the 
nonrotating sample; the small shift of absorption away from the Larmor region seen 
in Fig. IA3I arises from the local contribution Flh ■ 

In general, metastable states containing a cluster of any number of vortices N 
between zero and the equilibrium number iVoq can be observed in rotating "^He-B 
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(Fis. IAl|l . In such a case, the magnitude of the azimuthal counterfLow is 

/ ^ f 0, < r < i?c, , . 

where the cluster radius Rc = R^jN/N^ . The number of vortex Hnes in a given 
cluster can be determined by analyzing the counterflow peak height, either by 
comparing to numerically calculated equilibrium textures and their corresponding 
spectra |144| . or with purely experimental calibration techniques, which have been 
described in Refs. llHlj . However, in transient conditions, for instance after the 
turbulent burst, also more complicated flow profiles can arise, which are no longer 
necessarily purely azimuthal. This is the case with a twisted vortex cluster which 
is left behind by an advancing vortex front (Fig. I23|) . as discussed in Sec. [S] and 
Ref. |152| : the nonzero axial component of counterflow gives rise to excess absorption 
close to the Larmor frequency. Thus NMR in •^He-B can be used to distinguish 
between different flow states and vortex configurations, when good models for the 
corresponding counterflow proflles are available. Such models can be obtained from 
either analytical arguments, or from detailed numerical simulations of vortices. The 
latter ones are especially valuable in complicated dynamic situations. Thus, if the 
order parameter texture is known, the NMR absorption spectrum is obtained from 
Eq. ljA.8(l . However, the reverse is not generally true: a measured spectrum does not 
uniquely fix the order parameter texture, only once many more details are known 
about the measuring situation. 
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